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1. Introduction 

1.1. Objectives 

In this review article we focus on specific problems which are of interest both for 
numerical relativists and for geometers. A number of review articles has been devoted 
in the last years to the technical aspects and state of the art of each respective domain 
— e.g. Andersson [T^], Ashtekar & Krishnan Bartnik & Isenberg [53], Baumgarte 
& Shapiro [H], Berger [55], Booth [m], Cook [T57], Friedrich [530], Friedrich & 
Kendall Gourgoulhon 75711152], Gundlach & Martfn-Garcfa [MH], Klainerman 

& Nicolo |359j, Krishnan t370], Lehner [371] , Lehner & Reula [373], Rendall [i3S] . 
Reula [438], Shinkai & Yoneda [466], Thornburg [490], Winicour [515], and also [2]. 
More theme specific reviews will be referred to in the corresponding sections. The aim 
here is not the exhaustive description of the respective specific tools, but rather to 
identify and discuss conceptual or structural challenges in General Relativity which 
represent good candidates for a close collaboration between geometers and numerical 
relativists. 

This review is inspired in a workshop with a similar name that took place in 
the context of the General Relativity Trimester held at the Institut Henri Poincare, 
Paris, from November 20 to 24, 2006 |Jj, which brought together specialists from both 
the mathematical and numerical ends of General Relativity (GR). Due to the broad 
range of topics covered in the workshop we have taken the methodological decision 
to restrict our attention to numerical and mathematical aspects of the description 
of isolated systems in GR. This somehow reflects the underlying ultimate interests 
of the organisers, biased towards astrophysically motivated problems not including 
Cosmology. Unavoidably, not all the topics discussed at the time are covered in this 
review and on the other hand, some of the topics discussed in the review did not have 
a counterpart talk in the workshop. 

1.2. General thoughts 

1.2.1. Why from Geometry to Numerics? General Relativity is a theory that 
describes the gravitational interaction as a manifestation of the geometry of spacetime. 
It is therefore natural to expect that a geometrical perspective may be helpful, and 
even sometimes fundamental, for the further development of the theory at both 
conceptual and practical levels — the latter understood as the explicit construction 
of solutions. Does this point of view live up to the expectations raised? 

In our opinion, the answer to this last question has two facets. On the one hand, 
as discussed in |260) . understanding GR can be seen as tantamount to understanding 
the properties of its solutions, and in this sense "getting qualitative and quantitative 
(respectively analytical, theoretical, and numerical) control on the long time evolution 
of gravitational fields under general assumptions is still the most important open 
problem in classical general relativity". In this context there seems to be plenty of 
occasion for the successful interaction between geometry and numerics, not only in 
the direction suggested by the title of the review but also on the converse one, from 
numerics to geometry. 

Regarding the second facet, arguably, the other great challenge for the theory 
is to make full contact with observations — in particular, nowadays, in the area 
of gravitational waves physics. The nature of this specific endeavour requires the 
efficient calculation of spacetimes describing isolated systems. This calculation, to 
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be physically realistic has to be performed using numerical simulations of some 
particular formulation of the Einstein equations. In this context, the fruitfulness of 
using geometry as a guideline in numerical applications seems much harder to assess, 
since the understanding of a particular analytical aspect of a geometric problem does 
not ensure its successful numerical implementation. Furthermore, it should not be 
forgotten that one is talking about the (potential) interaction of two communities 
with a different history, traditions, languages and in great measure, objectives. 

1.2.2. Numerical Relativity and Relativistic Astrophysics. As somehow hinted in the 
previous paragraphs, it is natural to assume that a "generic" relativist would have 
interest in both studying the nature of solutions of theory — and hence improving his 
or her understanding of the theory — and also of making contact with the observation 
— which then would help to set the context and validity of the theory. Given the 
broad spectrum of research topics that relativists deal with, it is to be expected that 
the particular choice of research problem will set the emphasis on either conceptual 
or observational aspects. In particular, this applies to the case of a numerical 
relativist. For example, from a simplified point of view a "gravitational collapse 
code", can be potentially read as a tool to either understand conceptual issues like 
Cosmic Censorship or as a tool to predict astrophysical phenomena such as supernovae 
explosions. In the first case, it is natural for the numerical relativist to gain insight and 
even try to simplify the code through a geometric view point. In the second case, she 
or he has no option but to attempt to include all physical ingredients that are believed 
to play a role. Therefore, geometric insights may be comparatively less relevant. In 
any case, whether a geometry point of view is used or not reflects, ultimately, in the 
way the code is constructed and read. Summarising, a given problem in numerical 
relativity can be labelled as either geometric or as astrophysical numerical relativity. 
The tenet of this review is that the effects of the interaction between mathematical 
and numerical relativists is bound to be more fruitful in the former case, while in the 
latter the rewards would be more indirect. This duality will be a recurrent theme in 
the review. 

1.3. Tensions and synergy between Geometry and Numerics 

As a first contact with the topic of the review we present some general reflexions which 
illustrate the possible tensions and synergies to be expected as result of the interaction 
between geometry and numerics. 

Global versus local. Global issues dominate modern mathematics. On the other hand, 
most numerical simulations are local in time. This leads to a tension between the tools 
and goals in each community. An example of how different the points of view could be 
is given by the notion of black hole. The classical definition of a black hole makes use 
of global ingredients — e.g. [177] . On the contrary, everyday numerical simulations 
make use of quasi-local characterisations related to the notion of trapped region. This 
tension reaches its climax in statements like there are no results about the existence of 
asymptotically flat vacuum black hole spacetimes with radiation — see the discussion in 
12.21 — by a mathematical relativist that can provoke the smile of a numerical relativist. 
However, this tension is not as dramatic as it appears at first sight. On the one hand, 
it is very possible that numerical simulations of binary black holes could provide the 
missing clues for an eventual global existence proof of the existence of dynamical 



From Geometry to Numerics 



4 



black hole spacetimes |427| . On the other hand, the global framework provided by 
the Cosmic Censorship Conjecture acts precisely as guarantor of some of the technical 
tools used in numerical simulations — for example excision or punctures. 

Good geometric or analytic properties as a hope for good numerics. To a geometer 
who is inexperienced in the field of numerics, one of the first reality checks he or she 
has to confront is that analytic well-posedness and convergence of a numerical solution 
to the analytic solution, do not guarantee the long term stability of the simulation 
— this issue has been clearly discussed in ^S3]. This last point raised before the big 
break-throughs in the simulation of black hole binaries is still valid as regards to the 
realistic and accurate extraction of physics [4241 11411 159j . For this, it is possible that 
an analytic or geometry insight may provide the crucial ingredient. 

Geometry as a way of prescribing and extracting physics. Numerics relies on the 
use of coordinates. This results, unavoidably, in ambiguities in the extraction of 
physics if geometric notions are not employed. However, these methods are global 
and hard to implement numerically. In practical applications, one could use, for 
example, perturbative or post-Newtonian notions to extract the physics. Nevertheless, 
this approach already relies in the acceptance of a certain global behaviour of the 
spacetime. But more importantly, the recent advances in the numerical simulations of 
spacetimes could provide the opportunity for the reassessment of geometrical objects 
as tools for the extraction of physics. 

We conclude the section presenting some examples which give a taste of both the 
difficulties and rewards of the interaction between geometry and numerics: 

- Application of ideas from dynamical horizons to the extraction of physical 
parameters. The development of some quasi-local approaches to black horizons 
— see section 16.21 — has been directly inspired by the needs of the numerical 
implementations. In particular Dynamical Horizons have already provided first 
examples of successful interaction in the calculation of mass, angular momentum 
and associated parameters — cf. Krishnan's IHP talk and ^57 , ,454, i455t 1371} I146j . 

- Semiglobal evolution from small data. The existence results of the development 
of hyperboloidal initial data close to Minkowski by Friedrich |247j has been 
illustrated by the numerical simulations of Hiibner — see [3281 . 

- Hawking mass. The straight-forward implementation of some natural geometric 
notions meets often with numerical difficulties. This has been exemplified for the 
case of the Hawking mass in Schnetter's IHP talk [T]. 

- Calculation of high order derivatives. Any numerical simulation contains noise 
which tend to be amplified in the numerical evaluation of derivatives — see again 
Schnetter's IHP talk [Ij. This, for instance, complicates the implementation of 
constructions involving derivatives of the curvature tensor — like in the case of 
the metric equivalence problem in GR |355[ I356j . 

- Handling of coordinates. Sometimes the use of coordinates which are easy to 
handle from a numerical point of view, could lead to complications like coordinate 
singularities and/or singular equations — see e.g. Schnetter's IHP talk Ij. In some 
circumstances a detailed geometric analysis does suggest a reformulation of the 
problem which would solve this particular tension [28l [2^ . 
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- Penrose inequality. Classical bounds on energy loss in the head-on collision of 
black holes using the geometric Penrose inequality — see e.g. [276] — have been 
shown to be overtly optimistic by numerical simulations. On the other hand, lines 
of thoughts based on the Penrose inequality have been useful in the discussion of 
initial data [209] . 

- Geometry-inspired initial data sets. Attractive geometric properties have been 
employed and prescribed in the construction of initial data sets for black hole 
spacetimes — see among others e.g. |1961 11981 14371 1349] . However, the numerical 
usefulness of these geometric inputs is still to be assessed. 

- Geometry-inspired evolution systems. The close interaction between geometry 
and numerics is illustrated in the application of geometry as a guideline in the 
development and implementation of certain GR evolution systems — e.g. Z4 
formalism in |107j and also [44 1|, I444|, I523|. I524j . However, the role of Geometry 
is not so evident in most numerical codes — see section [SJ 

It may be of interest for the reader to look at the list of successes given in [12] . 
Gautionary notes 

Although this review stems from a workshop with a similar name |Tj , it should not be 
regarded as a proceedings or a systematic account of the former. We have used the 
diverse contributions to inspire a reflection on the interaction between geometric or 
analytic methods and numerical tools in the context of GR. In this sense, unavoidably, 
it projects our personal prejudices and it is certainly not comprehensive. We are 
indebted to all the participants for the effort put on their presentations and apologise 
beforehand in case of any misrepresentation — or omission — of their ideas. The 
presentations during the workshop will be referenced in a bundle as IHP talks. 

It is clear — although it is already evident from the introductory paragraphs — 
that the word Geometry has been used in a loose sense that encompasses more than 
Differential Geometry and includes, for example, (global and local) analysis, PDE 
theory or Group theory. 

As already mentioned, our discussion will be centred in numerical and 
mathematical aspects of isolated systems in GR. There are other streams of Relativity 
which have experienced a fruitful interaction between the analytical and numerical 
communities: namely, on the one hand higher dimensional spacetimes and on the 
other the study of Cosmological spacetimes. In the latter case, this interaction has 
been reported in the review [5S] — see also the discussion in [T^]. Without going into 
much detail, we just mention that the analysts (geometers) working in mathematical 
Cosmology have gained much from insights obtained through numerical experiments 
in which, for example, the behaviour of Gowdy and Bianchi IX spacetimes is explored 
towards the initial singularity — for a recent work on the subject see |93j . In particular, 
it was using this sort of numerical experiments that the existence of the so-called spikes 
— non-smooth behaviour and big gradients of certain field quantities as one approaches 
to the singularity — was firstly attested in the Gowdy spacetimes [88] . The analytical 
existence and nature of these spikes has since been thoroughly analysed. This is a 
prime example of the use of numeric techniques as a tool to explore the nonlinear 
behaviour of the gravitational field. 

We start from the premise that Mathematical Relativity and Numerical Relativity 
are research areas on their own right. This review is addressed to those researchers 



From Geometry to Numerics 



6 



working in any of these two communities who feel that there is place for a fruitful 
interaction. We assume a broad general knowledge of GR. We will not be extremely 
rigorous in our mathematical presentation and will refer the reader to the specialised 
literature. 

The article is structured as follows. Section [2] reviews some conceptual issues 
concerning the notion of isolated bodies in GR, including asymptotic flatness and 
the black holes. In section [3] the initial value problem in GR is presented as the 
appropriate approach to the generic construction of spacetimes. In particular, 13.11 
reviews the different types of initial value problem. Section [1] explores some aspects 
of the construction of solutions to the constraint equations which we believe are, or 
could be, of interests for numerical applications — the topology of initial hypersurfaces, 
the conformal method, the thin sandwich and the recent developments in gluing 
constructions. Section [5] discusses various aspects of evolutions formalisms. Here, 
particular emphasis is given to initial value problems of Cauchy type, although some 
observations concerning the characteristic and hyperboloidal problems are made. 
There is, also, a discussion of the conformal field equations. Section [6] reviews 
some aspects of the modern description of black holes including recent results about 
uniqueness and rigidity of the Kerr spacetime, quasi-local definitions of black holes, 
geometric inequalities, helical Killing vectors and aspects of "puncture" black hole 
evolutions. Section [7] discusses miscellaneous topics that did not find an adequate 
place in the structure of the review, but that we believe could not go unmentioned. 
Finally, [8] provides some conclusions and a brief list of topics where, we believe, the 
interaction between Geometry and Numerics will acquire particular relevance in the 
coming years. 

2. Isolated systems in General Relativity 

As mentioned before, this review will be concerned with solutions of General Relativity 
describing isolated systems. GR is formulated as a gauge theory of gravity, usually in 
terms of the Einstein field equations (EFE) 

Rfii^ - ^-R.9a"^ = SttT^^, , (1) 

where g^^ denotes the metric of the spacetime, R^^ its corresponding Ricci tensor, 
R the Ricci scalar and T^^ is the so-called stress-energy tensor — in geometric units 
c = G = 1 . Solutions to the EFE correspond to geometries rather than to metrics given 
in a specific coordinate system, since metrics related by spacetime diffeomorphisms 
are physically equivalent. As a consequence of this, equation ([T]) understood as a 
partial differential equations (PDF) system for the metric components g^i, constitutes 
an underdetermined (incomplete) system that cannot produce a unique solution. An 
extra ingredient, namely a gauge fixing procedure, must be provided to establish a 
standard PDF problem — note this remark applies to any geometrically introduced 
PDF. Here we adopt a broad notion of gauge which includes not only a coordinate 
choice but could also contemplate, for example, a tetrad choice and even perhaps a 
choice of conformal factor. 

We are interested in solving equation ([1]) subject to boundary conditions 
describing an isolated system for which the effects of the Cosmological expansion 
are neglected. A first way of approaching this problem is the construction of exact 
solutions to equation ([1]) by making assumptions on the symmetries of the spacetime 
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and/or on the algebraic structure of their curvature tensors — the most comprehensive 
treatise on exact solutions is |477j : reviews concerned with the physical aspects of 
exact solutions can be found in [571 HH] ■ This way of doing things has provided deep 
conceptual and physical insight, as it has supplied the tools of most of the observational 
predictions of GR. Furthermore, it is the basis of perturbative analysis. 

However, if one is interested in making statements about generic spacetimes 
having the right kind of boundary conditions and on how to construct them in a 
systematic manner, then one is led to consider an initial value problem or more 
generally and initial boundary value problem. Initial value problems can be of Cauchy 
type — that is, the initial data is prescribed on a Cauchy hypersurfac^ — or not — 
like in the case of the characteristic and the hyperboloidal [244j initial value problems 
— for a more complete discussion on these ideas see subsection 13.11 Only Cauchy 
initial value problems allow, in principle, to reconstruct the whole spacetime, while 
the characteristic and hyperboloidal problems allow at most to reconstruct the domain 
of influence of the initial hypersurface. 

The overall strategy is to prescribe suitable initial data on the initial hypersurface 
and then to evolve them by means of the EFE and the appropriate equations for the 
matter part. This is very much in the spirit of the approach taken by numerical 
relativists. There are, however, a number of caveats to this approach — like existence 
and uniqueness of solutions, stability and global issues including Cosmic Censorship — 
which are the concern of mathematical relativists (see |449| for a review addressing 
some open problems in mathematical relativity and articulated around the notion of 
Cauchy hypersurface). In particular, the satisfactory resolution of these caveats by the 
mathematical community is the guarantor of the consistency of the strategy employed 
by the numerical relativists. 

Very often, the needs of the numerical community are well ahead of the 
developments of the mathematical community. This applies, in particular, to the 
study of the motion of isolated bodies in GR, a problem of clear interest for 
astrophysicists which has been dealt numerically for a long time now — see e.g. 
|76| . The existence issue, that is taken for granted in the numerical community, is 
a challenging mathematical problem. We emphasise that this last point represents a 
crucial aspect in the overall coherence of the problem. The afore-discussed state of 
affairs is exemplified by the recent work on the existence of solutions describing the 
motion of isolated bodies — dust balls in this case — by Y. Choquet-Bruhat and H. 
Friedrich [161j . See also Choquet-Bruhat 's IHP talk [T. Finally it must be said that 
a converse phenomenon is also true. The numerics for many mathematical questions 
is still awaiting to be developed — for example the well-posed, constraint preserving 
initial boundary value problem constructed by Friedrich & Nagy |257j . 

2.1. Asymptotic flatness 

The physical intuition suggests that the spacetime describing an isolated body far 
away from it, should be Minkowskian in some appropriated sense. Such a behaviour 
is generically referred to as asymptotic flatness. As already mentioned, checks on the 
consistency of the theory requires a rigorous formulation of this physical intuition. 

t For a discussion of the relationship between global hyperbolic spacetimes and the existence and 
properties of smooth Cauchy surfaces, in particular the existence of Cauchy time functions, see 

[sniEiiiig. 
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The metric of an asymptotically flat spacetime, when expanded in terms of some 
convenient distance parameter, should coincide with the flat Minkowski metric at 
leading order. The physics of the problem is, on the other hand, encoded in the low 
order terms. In particular, one should expect to flnd there information about the mass 
and radiative properties of the system. Now, it is natural to ask, firstly, whether it 
is possible to rephrase the above statements in a coordinate independent way; and 
secondly, if one can make generic statements about the properties of the low terms in 
the asymptotic expansions of the metric of an isolated system. To this end, Penrose 
introduced the more specific notion of asymptotic simplicity — see e.g. |412[|4lS] for a 
precise definition. Roughly speaking, a spacetime {M.,g^i,) is said to be asymptotically 
simple if there exists a positive scalar $7 (the conformal factor) and a spacetime with 
boundary {M,g^^) — the compactified, unphysical spacetime — with g^i, = ^"^g^v 
such that the boundary of A^, defined by = 0, has a structure similar to that of the 
standard compactification of Minkowski spacetime — that is, the conformal boundary 
consists, at least, of a null hypersurface (null infinity) = ^~VJJ^'^ . In the particular 
case of the Minkowski spacetime, the conformal boundary consists, in addition, of 
three points: — spatial infinity — and i~ and i"*" — respectively, past and future 
timelike infinity. In general, these additional points will not be present in the conformal 
completion of a generic asymptotically simple spacetime. However, if the spacetime is 
constructed as the development of asymptotically Euclidean initial data — see equation 
([2|) — then «° will be part of the conformal boundary by construction. The situation 
with respect to i~ and is more complicated. If the relevant spacetimes contain 
black holes, then there is no reason to expect that be present in the conformal 
completion, or if they are, that their structure should be that of a point — this is 
a difficult open question which will benefit from input from numerical simulations. 
In any case, and to summarise, the definition of asymptotic simplicity gets around 
these issues by being only concerned with null infinity. The definition of asymptotic 
simplicity is geometric in what it avoids the use of coordinates. Penrose's original 
definition requires the boundary of the unphysical spacetime to be smooth (C°°). 
Penrose's original idea was to use the notion of asymptotic simplicity as a way of 
providing a characterisation of isolated systems in GR: the Jar fields of spacetimes 
describing isolated systems should be asymptotically simple in the sense that they admit 
a smooth conformal compactification to null infinity. The latter suggestion is known 
in the literature by the name of Penrose 's proposal. 

The conceptual and practical advantages of the use of the compactified picture 
to describe isolated systems have been discussed at length in the literature — see e.g. 
[36, 249, 250, 252, 253 J275j . Here we just mention that: i) the use of the compactified 
picture allows to rephrase questions concerning the asymptotic decay of fields by 
questions of their differentiability at the conformal boundary; ii) if the spacetime is 
asymptotically simple, then it is possible to deduce a certain asymptotic behaviour for 
the components of the Weyl tensor known as the peeling behaviour — see e.g. [236]. If 
the spacetime is such that it would not admit a smooth compactification — say, just 
C*^ instead of being C°° — then part of this formalism can still be recovered, but one 
has to be more careful — see e.g. [183) . 

The idea of asymptotic flatness and the notion of asymptotic simplicity are 
inspired by the analysis of exact solutions to the EFE. If one wants to establish more 
general properties of this class of spacetimes one has to resort to a formulation of the 
problem based on an initial value problem and then make use of analytic techniques to 
obtain qualitative information about the solution or construct the solution numerically. 
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Here, again, intuition suggests that the right initial conditions for obtaining an 
asymptotically flat spacetime are asymptotically Euclidean ones. More precisely, one 
could require an initial hypersurface E with at least one asymptotic end (a region 
which is topologically M'^ minus a ball) in which one can introduce coordinates such 
that 

where = ^J {x^)"^ + (a;^)^ + {x^Y , and 7y and Kij are, respectively, the first and 
second fundamental forms of S — see section |4] for more on this. More general 
notions of asymptotic Euclideanity have been considered in the literature — see e.g. 
|1521 11531 1173 L The one given here, has a well defined ADM mass at the asymptotic 
end in question. 

The first results concerning the existence of solutions to the vacuum EFE under 
the above boundary conditions can be traced back to [231] in which local existence 
in time was shown to hold — well-posedness. Choquet-Bruhat and Geroch introduced 
the notion of a maximal future development of an initial data set, and showed that for 
each given initial data set there is a unique maximal future development }162j . The 
work by Penrose has given a negative answer to the question of completeness in his 
singularity theorem [413j : given initial data where the initial Cauchy hypersurface is 
non-compact and complete, if the hypersurface contains a closed trapped surface — 
see section l6.2.2l — then the corresponding maximal future development is geodesically 
incomplete. A further milestone was the resolution of the so-called boost problem in 
GR given in [160|, I173j . Along these lines, the existence results of Friedrich |247j 
where it has been shown that the development of initial sufficiently small data given 
on an hyperboloid is complete. It should be noted that a hyperboloid is not a Cauchy 
hypersurface, and in this sense the results of \2A1\ are semiglobal. The analysis of 
the non-linear stability of the Minkowski spacetime of [172j has shown, among other 
things, the existence of vacuum radiative spacetimes which are asymptotically flat. 
More precisely, they show that every asymptotically flat initial data which is globally 
close to the trivial data gives rise to a solution which is a complete spacetime tending 
to the Minkoswki spacetime along any geodesic. There are no additional restrictions 
on the data. The result gives a precise description of the asymptotic behaviour at 
null infinity. Remarkably, however, the analysis is insufficient to guarantee the peeling 
behaviour of the solutions. An alternative proof has been provided by Lindblad & 
Rodnianski [375]. This provides a global stability result of Minkowski spacetime 
for the vacuum EFE in wave coordinate gauge — see section 15.11 The initial data 
coincides with the Schwarzschild solution in the neighbourhood of spacelike infinity. 
This result is less precise as far as the asymptotic behaviour is concerned. More 
recently, following the spirit of Christodoulou & Klainerman's proof, Bieri — see her 
IHP talk fTj — has considered more general asymptotically Euclidean initial data — i.e. 
with less decay and one less derivative — and shown the existence of a solution which 
is a geodesically complete spacetime, tending to the Minkowski spacetime at infinity 
along any geodesic. The solutions have finite ADM energy, but the angular momentum 
on the initial hypersurface is not necessarily well defined. The question answered by 
Bieri's analysis was what is the most general non-trivial asymptotically flat initial 
data set giving rise to a maximal development that is complete. The remaining open 
question in this programme is to find sharp criteria for non-trivial asymptotically flat 
initial data sets to give rise to complete maximal developments. 
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In what concerns the existence of spacetimes which satisfy the peehng behaviour, 
Klainerman & Nicolo [3601 1361] have shown that the peehng behaviour is satisfied by 
the development of a big class of asymptotically Euclidean spacetimes. This result is, 
however, not sharp for it can be seen that initial data for the Kerr spacetime — which 
satisfies the peeling behaviour as it is a stationary spacetime, see e.g. |212[ 1197] — is 
not contained in the class of initial data sets they considered. 

Regarding the related question of the existence of asymptotically simple 
spacetimes, it is now known that there is indeed a big class of spacetimes with this 
property — see [179j . The essential ingredients for the latter result are the semiglobal 
results of [247j , and a refined version of the initial data sets which can be constructed 
using the Corvino-Schoen method. This construction allows to glue fairly arbitrary 
asymptotically Euclidean initial data sets to a Schwarzschild asymptotic end — for 
more on this gluing construction see section 14.51 It has to be mentioned that the 
class of spacetimes obtained in |179j are very special because by construction, they 
do not contain any radiation near spatial infinity. The fact that the initial data 
is Schwarzschildean near infinity, implies that the Newman-Penrose constants of the 
spacetime vanish — the Newman-Penrose constants are absolutely conserved quantities 
defined in the cuts on null infinity |404|, [405]. This value propagates along the 
generators of null infinity, and it allows us to conclude that the Weyl tensor vanishes 
at future timelike infinity — see [2591 1256] . 

Given the state of affairs described in the last paragraph, the challenge is now 
to find the sharp criteria on an initial data set to obtain an asymptotically simple 
spacetime. Seminal work in this direction is |251] — see also [255] — where a convenient 
framework to discuss the behaviour of the region of spacetime where null and spatial 
infinity meet — the so-called cylinder at spatial infinity — has been introduced. It is 
important to point out that this work and related ones — see e.g. [498[ 14971 15001 1501] — 
are carried out in the conformally rescaled spacetime (unphysical spacetime) and 
make use of the so-called conformal field equations, a generalisation of the EFE which 
exploits the extra gauge freedom — a conformal one — that has been introduced via 
the notion of asymptotic simplicity. The programme started in [251] and continued 
in [2551 14981 14971 15001 1501] has showed that there is a hierarchy of obstructions to the 
smoothness of infinity which control the smoothness of the spacetimes at null infinity 
— see also Valiente Kroon's IHP talk fl\ . These obstructions are expressible in terms of 
parts of the initial data of the spacetime. A further analysis of the obstructions points 
towards the following conjecture: if the development of asymptotically Euclidean 
initial data admits a smooth conformal compactification (in the sense of Penrose) 
at both past and future null infinity, then the initial set is stationary near infinity 
— see [498] 1497] 1500] 1501] . If shown to be true, the latter conjecture would give a 
prominent role to the initial data sets which can be constructed by means of the 
Corvino-Schoen gluing construction. It is worth mentioning that the developments of 
conformally flat initial data sets like the Brill-Lindquist jl29| . Misner [400] . Bowen- 
York 1 124] and Brandt-Briigmann [125] , as well as the data used systematically in the 
numerical simulations of binary black hole mergers — see the end of section 16.4. H — 
will have a null infinity (if any) which is not smooth. However the differentiability of 
the conformal boundary seems to be enough to guarantee the peeling behaviour [502j . 

Given an arbitrary asymptotically Euclidean initial data set (S, 7^-, i^y), 
using the Corvino-Schoen construction one could construct another initial data set 
(E'j 7j'j-, K'^j) which coincides with the former inside a compact set, but which is exactly 
stationary near infinity. The development of (E, /lij, i^^) will possess a non-smooth 
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but possibly peeling null infinity. On the other hand, the development of (S', h'^^, K'^j) 
will be smooth. A natural question is the following: is there any substantial difference 
between the physics of the two developments? This is a complicated question which 
will only be answered in a satisfactory manner by a close collaboration between 
numerical and mathematical relativists. It is very likely that the latter issue will not 
have any relevance for numerical relativists interested in astrophysical applications. 
On the other hand, they provide a whole unexplored area for those who want to use 
the numerical methods to understand the nature of GR. 

2.2. Black holes 

A black hole is a region of no escape which does not extend out to infinity. Traditionally 
— see e.g. [310j — the definition is made using the conformal compactification of 
spacetime. Namely, if the spacetime admits a suitable null infinity and the causal past 
of null infinity is globally hyperbolic — i.e. it admits a Cauchy hypersurface — then 
one defines the black hole region like the spacetime minus the causal past of future null 
infinity. The boundary of the black hole is the event horizon. A spacetime containing 
a single black hole contains two natural boundaries, one is null infinity and the other is 
the event horizon — an inner boundary. The original picture arose from an analysis of 
spherical symmetry. According to the so-called weak Cosmic Censorship picture — cf. 
section [^751 — the end state of any generic asymptotically Euclidean initial data set will 
be a particular stationary black hole spacetime: the Kerr solution. Naked singularities 
can actually form, but the resulting spacetime are unstable |157[ 11701 1171j and non- 
generic. The rigorous definition of a black hole makes use of global ingredients: in 
particular it requires the existence of a complete J^+. As already mentioned, there are 
no global existence proofs of dynamical black hole vacuum spacetimes. The latter is a 
point deserving further clarification. There are results concerning the global existence 
of dynamical black holes in, for example, the spherically symmetric Einstein-scalar 
field theory — e.g. [193] for the formation of black holes from a regular past — or in 
the vacuum 5-dimensional setting. Regarding the latter example, see |194] where the 
orbital non-linear stability — i.e. with respect to perturbations respecting symmetry — 
of the Schwarzschild-Tangherini black hole is discussed, and thus the existence of a 
global, non-stationary black hole is established. Unfortunately, these examples are not 
relevant for astrophysics. The available examples of 4-dimensional vacuum dynamical 
black holes suffer from some pathology — like the Robinson- Trautman spacetimes 
discussed in jl74] with a complete future null infinity, but where past null infinity 
is incomplete — or are local — like in the discussion of "multi-black hole" spacetimes 
of [184] ■ or that of spacetimes with isolated horizons of [375) . 

Among other things, a proof of the nonlinear stability of the Kerr spacetime would 
show the existence of a wide class of spacetimes with non-stationary black holes. The 
idea behind such a proof is to show that initial data which are in some sense close to 
Kerr initial data will give rise to a black hole spacetime whose global structure will be 
similar to that of the Kerr spacetime. In addition to proving the existence of dynamical 
black hole spacetimes, such a result would be an important step towards a proof of 
the (strong) Cosmic Censorship Conjecture — see section [3l A particular situation 
in which the latter conjecture has been proved in an asymptotically flat context is 
the spherically symmetric Einstein-scalar field system studied by Christodoulou — 
see |1661 11671 11681 11691 I171| . A crucial property of stationary black holes is the 
so-called rigidity: that is, the fact that stationarity together with a couple of further 
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assumptions on the behaviour of the horizon imply further symmetries in the spacetime 
— i.e. axial symmetry. The ideas behind the rigidity theorems have been recently 
revisited by Isenberg & Moncrief with the aim of obtaining generalisations which are 
valid for higher dimensions — cf. Isenberg's IHP talk ^Ij and section IB. 1.1 1 

2. 3. Critique and practicalities 

There has been a number of critiques to the notion of asymptotic flatness, most 
notoriously the ones articulated by Ellis [2221 [223J . This critique has Cosmological 
motivations and arose from the desire of analysing the nature of the interaction 
between local physics and boundary conditions in the expanding Universe. In 
particular, there was a desire to understand what difference does the use of particular 
boundary conditions make on results depending in an essential manner on asymptotic 
flatness, like the peeling behaviour, the positivity of mass — cf. section l673l — and black 
hole uniqueness. Ellis argues, for example, that it is unrealistic when considering local 
gravitational collapse to worry about observations with infinite life-times. The latter 
is what is implied in the usual definition of a black hole — see however section 16.21 
The alternative to null infinity as a natural boundary of spacetime given in [222j is to 
consider a world-tube — finite infinity — with a timelike boundary located at a spatial 
radius which is sufhciently far away from the sources. If brought to the context of 
an initial value problem, Ellis proposal implies the use of an initial boundary value 
problem to describe the physics of isolated systems in which the effects of the Universe 
are put in by appropriate boundary conditions. The well possedness of the initial 
boundary value problem has been firstly discussed in |257j . The analysis given there 
shows that if one makes use of a finite infinity then there is no covariant way of 
prescribing the boundary information. This, in turn, would lead to ambiguities in the 
extraction of physics. One of the rationales behind the introduction of the null infinity 
formalism was precisely to avoid this sort of problems. For a neat discussion of these 
issues, with an emphasis to the extraction of radiation see [234j . 

This critique is remarkable in what that if one takes it away from its Cosmological 
context it could also be employed when addressing the relevance of global issues for 
numerical relativity. Most numerical simulations make an implicit use of the notion of 
finite infinity by calculating numerically the solution to an initial boundary problem. 
In this approach there is always the potential of "messing up" the physics of the 
problem by setting boundary conditions which are not appropriate — see section [5AT] 
There are, however, also simulations which compactify in space, see [427) . Moreover, 
the numerical simulations do not cover an infinite time although given the current 
advances in the field there seems to be no essential difficulty for letting the simulations 
"to run as long as it is necessary" . 

It is worth noticing that there have been some semiglobal numerical calculations 
using the conformal field equations which exemplify the hyperboloidal existence results 
of [HHI —see [323 1323 WM WM- There are some further attempts to calculate 
portions of the Schwarzschild from hyperboloidal data — [330] unpublished. For a 
critical review of the successes and problems of this programme see [332] . More 
recently, there have been some global calculations of the Schwarzschild spacetime 
[523[ I524j . These calculations make use of a more general version of the conformal field 
equations and conformal Gauss coordinates. These numerical simulations exemplify 
the analysis of this type of coordinates given in [254] . We also mention the recent 
fully pseudo-spectral scheme [318| for solving hyperbolic equations on conformally 
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compactified spacetimes. 

The discussion in the previous paragraphs leads to a situation where besides 
a number of simulations tailor-made for exemplifying some analytical aspects of 
solutions of the EFE, a majority of the state-of-the-art numerical codes used to 
simulate spacetimes containing dynamical black holes are not designed for addressing 
global issues in GR. So, what is the relevance for the numerical community of the 
work on global issues carried out by mathematical relativists? As mentioned in the 
introductory section, the role is to ensure that some of the approaches implemented 
in numerical simulations are justified. For example, if the non-linear stability of the 
Schwarzschild turned out to be false, then it could happen that the development of 
an asymptotically Euclidean black hole initial data would not have the asymptotic 
structure of the Kerr spacetime, and thus measurements of gravitational radiation 
at finite radius would be rethought. Similarly, if the no- hair theorems were not 
valid, there would be no justification for employing perturbative methods on a Kerr 
background to analyse the late stages of the evolution of a black hole binary merger. 
These examples are extreme but highlight the need of having the theoretical framework 
on a sound ground. 

Finally, we would like to retake the point raised in [12] that the problem of the 
stability of the Kerr spacetime opens a natural arena for the interaction of numerical 
and mathematical relativity. Indeed, it is easy to imagine that tailor-made simulations 
could provide crucial insights concerning the asymptotic decay of the gravitational field 
near the horizon. This interaction would be, notwithstanding, on the lines of from 
Numerics to Geometry. 

3. Initial value problems 

If one wants to discuss generic properties of solutions to the Einstein field equations 
which describe in some appropriate sense isolated systems, one has to find a 
systematic way of obtaining these solutions. The only approach to the construction of 
solutions which seems general enough is that of an initial value problem (IVP). This 
mathematically sound approach is perhaps not ideal from a physical point of view, 
and indeed a number of objections have been raised. One could argue that all major 
observational predictions of GR have been obtained by means of the analysis of exact 
solutions. 

A first (conceptual) objection consists in making sure that one is able to 
reconstruct, starting from initial data, the whole maximal extension of the spacetime 
— the strong Cosmic Censorship |4161 14171 [TT| . This has implications on causality 
violations and predictability of the theory — see e.g. [449]. 

Another objection questions to what extent an approach to GR based on an 
initial value problem is appropriate to obtain predictions which can be contrasted 
with measurements. One of the crucial problems behind is that for a certain system of 
physical interest, there may be many possible ways of constructing physically plausible 
initial data. This is particularly the case when attempting to construct initial data for 
spacetimes containing dynamical black holes. In order to construct data for a black 
hole spacetime, one should be able to find a way of identifying the black hole in the 
initial data. For these purposes, the notion of event horizon is of no use for it is a 
global quantity which is only known once the whole spacetime — up to null infinity — is 
known. Instead, one can use the notion of apparent horizon: for if there is an apparent 
horizon on the data it must be contained in the black hole region — see section [621 In 
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other words, the presence of an apparent horizon indicates the existence of a black hole 
— if the weak Cosmic Censorship is true. Furthermore, data with an apparent horizon 
will be geodesically incomplete as a consequence of Penrose's singularity theorem [413] 
— thus, indicating the formation of a singularity. It should be emphasised that a black 
hole data does not need to contain an apparent horizon, but if it contains one, it must 
be a black hole. Similarly, it is not clear whether data with two apparent horizons will 
contain two black holes. However, this is a reasonable assumption to begin with. To 
make things worse, note that many of the approaches for constructing initial data for 
black holes are explained more in terms of mathematical advantages, like the simplicity 
of the equations — like in the case of conformal flatness — than in physically based 
considerations — see however section 16.21 for an account of the ideas introduced by 
the quasi-local approach to black holes. To summarise, physics seems to be hard to 
encode in initial data sets for the EFE. 

There has been a number of works devoted to compare the results of numerical 
simulations in the search of robust aspects — that is, structures in the numerically 
calculated spacetimes which are independent of the particular type of initial data 
being used — see e.g. [60l 14741 149 Ij and IHP talks by Hannam and O'Murchadha. 
Crucially, the usefulness of numerical simulations as a source of wave form templates 
for the detection of gravitational radiation in the current generation of interferometric 
gravitational wave observatories depends on the assumption of the robustness. From 
the point of view of a mathematical relativist or a geometer, this issue is essentially 
open and unexplored. The problem is certainly complicated, and possibly one will have 
to wait several years before some qualitative statement can be made. The resolution of 
the non-linear stability of the Kerr spacetime may provide first insights and rigorous 
answers to the question of robustness of initial data sets. 

3.1. On the different types of initial value problems 

Initial value problems can be classified according to the nature of the hypersurface 
on which the initial data is prescribed. Usually, the initial hypersurface E is taken 
to be spacelike. As mentioned in section 12. H when discussing initial value problems 
for isolated systems the spacelike initial hypersurface S is naturally assumed to be 
asymptotically Euclidean. This particular type of initial value problem in GR will be 
referred to as a Cauchy initial value problem as it would allow, in principle, to recover 
the whole spacetime — see references |3101I504] for a rigorous discussion of these issues, 
and ^258J for a recent review on the subject. The seminal work of Choquet-Bruhat [231] 
has shown the well-posedness of the Cauchy initial value problem in GR: an existence 
theorem local in time — that is, existence is guaranteed for a small time interval, in 
a neighbourhood of the initial hypersurface. The portion of spacetime recovered by 
this local existence theorem is called the development of the initial data. At this level 
is not possible to discuss global existence of solutions without more assumptions or 
information on the initial data — for further discussion on this, see for example the 
review |359) . Nevertheless, it has been shown that the development has a unique 
maximal extension [162] . 

Initial value problems on spacelike hypersurfaces are not necessarily Cauchy initial 
problems. For example, the hyperboloidal initial value problem first introduced by H. 
Friedrich in [244j . on which initial data are prescribed on a spacelike hypersurface 
— which could be thought of as intersecting null infinity — is not of Cauchy type. An 
initial value problem on a hyperboloidal hypersurface would not allow us to recover 
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the whole spacetime. Moreover, it is not possible to find out where the development 
of hyperboloidal data lies in a globally hyperbolic spacetime, more specifically, one 
cannot relate the hyperboloidal data with a hypothetical Cauchy data. 

As pointed out in [478] . the Cauchy problem in GR is in some sense natural for 
theoretical discussions, but strictly speaking it does not necessarily correspond to the 
way things are done in numerical relativity, where one usually aims to solve the EFE 
only within a compact domain — cf. sections 12.31 and 15.4.11 The latter requires the 
solution of an initial boundary value problem^ where in addition to the data to be 
prescribed on a spacelike domain, one also has to prescribe extra data on a timelike 
boundary. The presence of this extra data — the boundary data — raises the question 
of which data can and must be specified on the timelike boundary in order to have 
a solution — and further, whether there are any compatibility conditions required 
between the initial data and the boundary data. A first rigorous discussion of the 
boundary value problem in GR, including a local existence result has been provided 
by Friedrich & Nagy |257| . A recent IVBP well-posedness result in the context of the 
harmonic formulation — see section 15. H — has been presented in [368[ I367j . 

The initial data must not necessarily be prescribed on spacelike hypersurfaces. 
One could also prescribe it on null hypersurfaces, in which case we speak of a 
characteristic initial value problem. However, as the domain of influence of a non- 
singular null hypersurface — i.e. a null hypersurface without caustics — is empty, one 
needs, for example, to consider data prescribed on two intersecting null hypersurfaces 
in order to obtain a non-trivial development. A local existence theorem for this sort 
of problem has been given in [435j . A variant of the problem is to prescribe data on 
a light cone. There is no existence result for the Einstein equations in this case — see 
however [246] . 

The characteristic initial value problem in GR can be traced back to the seminal 
work of Bondi and collaborators on gravitational radiation in the early 1960's — 
see |116l 1447] . An early systematic discussion of the characteristic problem can be 
found in [448) . From the point of view of asymptotics, one can formulate the so- 
called asymptotic characteristic initial value problem in which initial data is prescribed 
on both an outgoing light cone and on null infinity. The portion of spacetime to 
be recovered from this initial problem lies at the past of these null hypersurfaces. 
Existence results for this type of characteristic initial value problem have been given 
in |242|, 12411 1243( I350j in the context of the conformal field equations. An analysis 
of the well-posedness of the problem more in the lines of the original work of Bondi 
and collaborators can be found in [262j . Another type of characteristic initial value 
problem is the 2-1-2 initial value problem — see e.g. |217j — and the initial value 
problems with data prescribed at past null infinity of [2461 1248j . 

On a characteristic initial value problem, a part of the EFE reduce to interior 
(transport) equations on the null hypersurface which can be recast as ordinary 
differential equations along the generators of the null hypersurface. Further, the 
data are free — i.e. not subject to elliptic constraints like in the case of spacelike 
hypersurfaces, see section |4l This feature was crucially exploited in the seminal work 
on gravitational radiation. The characteristic problem is naturally adapted to the 
discussion of gravitational radiation, however, it tends to run into problems when 
caustics form. Further, physics is hard to prescribe on null hypersurfaces. 

Finally it is noted that an IVP of a mixed type on which data is prescribed on 
a combination of spacelike and null hypersurfaces has been considered: the Cauchy- 
characteristic or matching problem — see [94] and the review |515| . 
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3.2. Gauge reduction 

As it has been discussed at the beginning of section [21 the formulation of the EFE 
as a standard PDE problem requires the adoption of a gauge choice. This reduction 
process — see e.g. |260[ 1258) — involves four different differential systems: 

i) The main evolution — also, reduced or relaxed — system, whose solution provides 
the spacetime geometry in a certain coordinate system. 

ii) The system of constraints, which controls the permanence in the submanifold of 
solutions of the theory. 

iii) The gauge system, which allows the fixation of a coordinate system and therefore 
casts the main evolution system as a genuine PDE system. 

iv) The subsidiary system, consisting in the evolution equations of the auxiliary 
(subsidiary) constraints introduced in the previous steps and guaranteeing the 
overall consistency. 

The different manners of dealing with the previous points define a variety of evolution 
formalisms in GR, where relative stress to be placed on systems i)-iv) depends on the 
specific aspect of the EFE one wants to address. In section [5] we will briefly review 
some of the more relevant evolution formalisms considered in numerical relativity with 
a focus on formulations based in the GR Cauchy IVP. Prior to this, one must discuss 
the system of constraints on a spacelike hypersurface. 

4. Initial data for the GR Cauchy problem 

The initial data to be prescribed in this IVP problem are given by a pair of symmetric 
tensor fields (7^ , i^ij) on the initial E, respectively corresponding to the induced 
Riemannian metric and the extrinsic curvature under the embedding of S into 
(M, g^j/). Such an interpretation of S as a spacelike hypersurface of (M, g^^^) demands 
the fulfillment of the so-called Hamiltonian and momentum constraints, respectively 



associated with 7^^ , K is the trace of Kij , p is the energy density and J' the current 
vector. An initial data set is said to be maximal if ii' = 0. If in addition Kij = 0, 
then it is called time symmetric. 

We proceed now to discuss some aspects concerning the way the Einstein 
constraints ^ and ^ are solved, putting some emphasis on the relation of these ideas 
to numerical implementations — a major review on general aspects of the constraint 
equations can be found in ^65j; see also |,187ll420ll289j for a numerical counterpart. For 
concreteness we shall concentrate our discussion on the vacuum constraint equations. 

^.1. Topology of the initial hypersurface 

As mentioned already in section [21 from a Cauchy initial value problem point of view, 
the natural requirement to obtain a spacetime describing an isolated body is to have 
an initial hypersurface S that has at least one asymptotically Euclidean end — see 
equation The theory of hyperboloidal hypersurfaces is less developed, although 
some of the essential ideas of the asymptotically Euclidean setting can be retaken — see 
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again [551 HSl 1131 HH HI] • Initial hypersurfaces with more than one asymptotic end are 
said to have a non-trivial topology. Note however, that the introduction of multiple 
asymptotic regions is not the only way of having initial data sets with non-trivial 
topology. The introduction of a wormhole connecting two points in a hypersurface 
would also result in a non-trivial topology. This was firstly exemplified in [399j ; with 
the use of the Isenberg-Mazzeo-Pollack (IMP) gluing techniques — see section IT4l — it 
is now possible to construct this type of initial data in a systematic way. 

Initial hypersurfaces with non-trivial hypersurfaces have been routinely used in 
the mathematical literature as a way of manufacturing initial data sets for black 
hole collisions. The simplest example of the latter are the usual time symmetric 
initial data for the Schwarzschild spacetimes in isotropic coordinates, which contain 
two asymptotic regions connected by a so-called Einstein- Rosen bridge, throat or 
wormhole. This example can be readily be generalised to contain n -\- 1 asymptotic 
regions rendering the so-called Brill- Lindquist initial data |129j . Another approach, 
which renders time symmetric initial data on a hypersurface E containing 2 asymptotic 
regions and n wormholes and can be interpreted as the data for n black holes — the 
Misner initial data — was given in [400i . 

A first step towards the construction of non-time symmetric data for black 
holes was given by Bowen & York in |124j . where an analytic solution of the 
momentum constraint in the conformally flat case was presented, the Bowen-York 
second fundamental form — a discussion of the general solution of the so-called 
Euclidean momentum constraint can be found in e.g. |2051 l85j. The so-called Brandt- 
Brilgmann puncture initial data |125j can be in some ways regarded as the non-time 
symmetric generalisation of the Brill-Lindquist data. It is important to mention that 
all the initial data sets discussed in this last paragraph are conformally flat — see 
section [4.21 But this is merely a simplifying assumption and in no way essential. 

The fall-off conditions for 7^ and Kij to be asymptotically Euclidean have 
already been discussed in section 12.11 — see equation ([2]). As in the discussion of 
boundary conditions for the whole spacetime, the discussion of the fall-off conditions 
on E can be geometrically reformulated in terms of a conformally compactified 3- 
dimensional manifold E — see [551 IMl [2461 1251] . The compactification procedure is, 
in this circumstance, analogous to that induced by the introduction of stereographic 
coordinates to compactify M'^ to render the 2-sphere From this point of view, the 
compactified initial hypersurface E contains, say, m singled out points {zi, . . . ,im\ 
corresponding to the m asymptotic ends of the physical initial hypersurface. For 
example, the standard Schwarzschild initial data render a E with the topology of the 
3-sphere §^ and two singled out points. The Brill-Lindquist data for 2 black holes 
would render a S with 3 singled out points, while in case of the Misner data the 
resulting E has the topology of a torus with 2 singled out points. 

From a mathematical point of view, working with E instead of E has several 
technical advantages. It is simpler to prove existence of solutions for an elliptic 
equation on a compact manifold than on a non-compact one. The price one has 
to pay is that the equations will be singular at the points at infinity 
However, these singularities are mild. It is also simpler to analyse the fields in 
terms of local differentiability in a neighbourhood of the points at infinity than in 
terms of fall-off expansions at infinity in E. This approach to discussing the theory 
of the constraint equations has been introduced and used in various applications in 
[2461 [2511 [T96l[T97l[T98l [82]. 

The puncture initial data used in many of the simulations of the coalescence of 
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black holes, as introduced in [125] , is a sort of compromise between using the physical 
hypersurface E and the compactified one S. For punctured data, one of the asymptotic 
ends — the so-called reference, or physical asymptotic end — is not compactified. The 
resulting initial hypersurface S will have the topology of S"^ \i i . The singular behaviour 
of the constraints near the punctures reflects in a singular behaviour of their solutions 
that is mild and well understood — it is related to the singular behaviour of the Green 
functions of the elliptic operators appearing in the constraint equations — see [83l [84] . 
Numerical relativists interpret this phenomenon as a coordinate singularity at the 
puncture — see [531 [HI 13051 1306j . The attractive of puncture initial data sets for 
numerical relativists is that they represent black holes in R'^ without excision and it 
is well understood how to construct puncture initial data for any number of boosted, 
spinning black holes |125|. I199j . 

The relevance for numerical relativists of initial data sets with non-trivial topology 
stems from the fact that they allow to model black hole spacetimes by only making 
use of the vacuum EFE: indeed, if there is a wormhole in the initial data, then it can 
be concluded that the spacetime will be geodesically incomplete |265| . Note however, 
that in all this argument one is indirectly assuming that Cosmic Censorship holds and 
that the resulting spacetime is really a black hole spacetime. In particular, from a 
physical point of view one is not interested in what happens in the extra asymptotic 
ends — see however the discussions in section [B.4.2l on why the puncture methods work. 
From the point of view of a physicist, it is fair to say that the use of initial data sets 
with non-trivial topology is a trick — there have been, however, some works on trying 
to analyse the possible physical differences between initial data sets with different 
topology but describing the same system —see e.g. |283| I284[ I496| 146711 147j . 

4-2. The Conformal Method. 

A standard approach to study the space of solutions of the constraint equations is the 
so-called conformal method — or conformal Ansatz. The method not only provides 
a systematic approach to make statements about the existence of solutions to the 
constraint equations, but also gives a starting point for the numerical computation of 
the solutions to the equations. A recent account of what is known about solutions to 
the constraint equations is given in [65] — see also |187[ I420[ [28 9 j . 

The conformal method is based on the conformal rescaling of the metric 

7y = *^7y , (5) 

and on the combined conformal rescaling and tensor splitting of the traceless part 
A^^ of the extrinsic curvature K"^^ into a transverse-traceless and longitudinal part. 
Different — non-equivalent — constructions follow from the non-commutativity of the 
conformal rescaling and tensor splitting process. A key point in all of them is the 
need of a specific rescaling for the (conformal) traceless part of the extrinsic curvature 
in order to set the constraint equations in a form which is optimal for their analytic 
study. In particular, if K is constant, the conformal method reduces the problem of 
solving the constraint equations to a decoupled underdetermined elliptic problem. 

Much has been learned about solutions of the Einstein constraint equations using 
the conformal method. However, many other interesting theoretical questions cannot 
be addresses using these methods. The main problem being that the conformal 
method produces solutions of the constraint equations, essentially from scratch. It 
turns out that it is very difficult to encode physics into the conformal metric jij 
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and in the symmetric transverse-traceless tensor A^^ — which act as seeds of the 
method. For example, there seems to be no way of constructing systematically initial 
data for stationary solutions, although there is an a posteriori check for the case 
of time symmetric data [200] . Most of the natural simplifying assumptions that 
are directly suggested by the structure of the conformal method happen to have no 
essential physical or geometrical content — most notably the assumption of conformal 
flatness. This is bluntly illustrated by results showing that there are no conformally 
flat hypersurfaces in the Kerr spacetime [267| 14991 1498j . Some examples of solutions 
to the Einstein constraints where it has been possible to put some desirable physical 
properties with some level of success can be found in for example [196[ 11971 1209] . In 
these last references either a non-flat metric or a different choice of the Bowen-York 
Ansatz for the extrinsic curvature — see |124j — have been used. Another approach 
to prescribe physics is the use of tensors 7ij and Aij motivated by post-Newtonian 
analysis as in [408j . 

It is noted that given the relation of Kij with the canonical momentum in the 
Hamiltonian formulation of the theory [35] 

= ^{Kf^ - K'^) , (6) 

this view of the constraints can be referred to as a Hamiltonian point of view |520il421] . 
No reference to the slicing is present and in particular no reference to the lapse function 
or shift vector emerges in the formulation of the constraints. This formulation reflects 
the character of the initial data problem as a one-hypersurface embedding problem. 
In this context, in particular in the numerical community, the conformal Ansatz is 
also referred to as the extrinsic curvature approach. 

4-3. The conformal thin sandwich method. 

The idea of a thin-sandwich approach to GR — in the spirit of the Jacobi principle in 
Mechanics — seems to have been firstly introduced by Baierlein, Sharp and Wheeler 
(BSW) [5S]. They conjectured the possibihty of specifying the spacetime metric g^,ij 
on two hypersurfaces and then of recovering the spacetime in between. In particular, 
when applied to two infinitesimally closed spatial slices — the BSW Thin Sandwich 
(TS) conjecture — this approach would imply a manner of providing initial data for 
a Cauchy evolution of the Einstein system. The TS conjecture, through its two- 
hypersurfaces Ansatz, proposes a different perspective to the constraint problem as 
compared to the one-hypersurface embedding of the extrinsic curvature approach. In 
this setting one aims at prescribing the induced metric and its time derivative (7^ , 7ij), 
in what can be seen as a Lagrangian formulation of the problem. However, although 
this conjecture can be seen to work for particular examples, the analysis in [67 shows 
that the conjecture fails in generic cases. 

More recently — see also previous works |3411 13401 15141 1189j and section [B.4.H — a 
new interpretation of the Einstein constraint equations, partly on the spirit of the TS 
view point, has been given in |520j . In this approach, the freely specifiable data on the 
initial hypersurface S consists of a conformal metric jij and its time variation Uij = 7^ 
— which is imposed to be traceless. The latter can be shown to amount to the choice 
of the conformal metric in two infinitesimally close slices. Using the kinematic relation 
between the time derivative of the physical 3-metric and the extrinsic curvature, the 
traceless part Aij can be written explicitly in terms of the lapse function N and 
shift vector /3* — cf. notation in section 15.21 A function N conformally related to 
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the physical lapse is identified as the naturally freely specifiable parameter. This 
conformal rescaling together with the absence of rescaling of the shift vector, permits 
to derive the key conformal rescaling for Aij introduced by Lichnerowicz 1376] and 
leading to the decoupling of Hamiltonian and momentum constraints in the constant 
K case. The specification of the free data (7^^ , 7^ ; -fC, A^) transforms the constraints 
into an elliptic system on 5" and similar in form to the one obtained in the extrinsic 
curvature approach [520| . In this Conformal TS approach (CTS) the role of every part 
of the metric is explicit. The formulation has been extended and improved in [421j . 
which in particular reconciles the Hamiltonian and Lagrangian points of view. The new 
key feature is the introduction of a weight function a which permits to recover the 
commutativity between conformal rescalings and tensor splittings. Moreover, with 
the choice a — 2N, the vanishing of the tranverse-traceless part of the extrinsic 
curvature Al^j, characterises stationarity, in agreement with the general (intuitive) 
understanding of Al^j, as the piece encoding the radiative — or dynamical — degrees of 
freedom. Regarding the free data, the approach in [421] pushes forward the Lagrangian 
analogy. Understanding K as a configuration variable, it is proposed the prescription 
of K instead of N at the (first) initial slice. This translates into a fifth elliptic equation 
for the conformal lapse N which together with the CTS equations gives rise to the 
Extended Conformal Thin Sandwich (XCTS) system — see e.g. [4211 14201 170] for the 
list of the equations. 

One of the advantages of the CTS approach is that using it, it seems easier 
to prescribe physics than in the standard approach based on the extrinsic curvature 
approach. In particular, as discussed in jl871l289] . the CTS method motivates definite 
Ansatze for the free part in the construction of stationary initial data — at least in the 
axisymmetric case. More generally, this CTS approach has been used also to construct 
binary initial data sets in quasi- equilibrium — see e.g. [691 1114[ 12911 1495] or Uryii's 
IHP talk 1 for the case of neutron stars; [2901 MB ^Hl [US [HI] for the case of 
binary black holes; and [72l 14871 1294j for the mixed black hole-neutron star binary 
systems. 

4-3.1. Uniqueness issues in the XCTS system. In contrast to the CTS and closely 
related elliptic systems — see e.g. [4101 14111 I518| and also further references in 
[5201 142 1| — there is to date a lack of systematic analyses of the mathematical 
properties of XCTS system. In particular, no definite results are available about 
existence conditions in this system, and very little is known about uniqueness. 
Numerical evidence of the non-uniqueness of solutions of the XCTS elliptic system 
has been found [422] . This evidence is of relevance as the XCTS system or somewhat 
related ideas are used to construct initial data for binary neutron star and black hole 
binary systems — cf. the literature listed in the previous paragraph. Moreover, it 
also has implications for the groups implementing constrained evolution schemes — cf. 
section 15.2.21 

A first analysis of the non-uniqueness issue, looking in particular at the question 
of genericity, has been given in [7Dj . The most complete analysis available of the non- 
uniqueness of the solutions to the XCTS system, using Lyapunov- Schmidt bifurcation 
theory, has been given in [507| — see also Walsh's IHP talk [l]. It is found that if the 
linearised system has a Kernel for sufficiently large data, then there exists a parabolic 
branching of the solutions as the one happening in [422]. The prosaic reason for 
the existence of the non-trivial Kernel is the presence of a wrong sign in one of the 



From Geometry to Numerics 



21 



equations of the XCTS system — more precisely, in the equation for the lapse. Thus, 
standard methods based on the use of a maximum principle cannot be employed. 
The numerical evidence shows that the assumption of a nontrivial Kernel is certainly 
not far-fetched, and actually occurs in numerical implementations. The situation gets 
even worse if a treatment of boundary conditions is required — see e.g. |345j . 

The issue of the non-uniqueness of solutions to the XCTS equations is an example 
of the tension underlying the relation between numeric and mathematical relativists: 
a successful method of prescribing physics into numerical simulations with non- 
desired mathematical properties. In the absence of better alternatives, the numerical 
community seems prepared to take the risk. 

4- 4- Construction of initial data sets using the Isenberg-Mazzeo- Pollack gluing. 

In recent years it has been possible to address a number of issues concerning properties 
of initial data sets using a tool of geometric analysis: gluing — cf. Pollack's IHP talk. 
While the conformal method is a procedure to construct solutions to the constraint 
equations from scratch, gluing constructions allow to obtain new solutions by means 
of direct sums of two already known solutions. Analytic gluing techniques have played 
a prominent role in many areas of differential geometry: some notable applications 
include the study of the topology of 4-manifolds and the study of minimal and constant 
mean curvature surfaces in Euclidean 3-space. 

In [3371 13381 1182] a particular gluing construction for solutions to the Einstein 
constraint equations was developed under the assumption of a constant K — the 
Isenberg-Mazzeo-Pollack (IMP) gluing. This construction allows to demonstrate how 
spacetimes can be joined by means of a geometric connected sum, or how a wormhole 
can be added between two points in a given spacetime — at the level of the initial 
data. This construction makes it possible to address a number of issues concerning 
the relation of the spatial topology to the geometry of solutions of the constraints and 
the constructibility of multi-black hole solutions jl84j . 

The sharpest possible gluing theorem for the Einstein constraint equations has 
been developed in [182j . Using it, it is possible to show that for a generic solution 
of the constraint equations and any pair of points in this solution, one can add a 
wormhole connecting these points to the solution with no change in the data away 
from a neighbourhood of each of these points. Further, one can show that for almost 
any pair of initial data sets — including, say, a pair of black hole data sets, or a 
cosmological data set paired with a set of black hole data — one can construct a new 
set which joins them. 

In order to assess the relevance of the IMP gluing construction for numerical 
relativity, we present a rough overview on how the IMP gluing method works. Consider 
two initial data sets (Ei, 7^?^ , if -^) and {Ti2,"/fj, Kfj) for the EFE. Let pi € Ei and 
P2 S S2 be two arbitrary points. The idea of the IMP construction is to find an initial 
data set {1^12, jfj, KlJ) such that E12 has the topology of Si#I]2, where # denotes the 
connected sum of manifolds^ The pair K]^) is such that on Si \ {ball around pi\ 
it coincides with {jjj,Klj) and on E2 \ {ball around P2} with {jfj,Kfj). 

In order to construct the pair {'-flj,Kl^) satisfying the Einstein constraints, one 
considers a conformal metric jf^ with the property that on Ei \ B{pi) and E2 \ ^(^2) 

X The connected sum is performed by removing B{pi) = {ball around pi} and B{p2) = 
{ball around P2} and then joining Si \ B{pi) and 52 \ B{p2) at dB{pi) and dB(p2) with a cylinder, 
X, of topology 9B(pi) X [0, 1]. 
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coincides with and , respectively, and on the cyhnder T it interpolates between 
them. Similarly, one constructs a ^IJ-tiaceXess tensor Al'j which again interpolates 
between A\j and A^y By construction, the resulting pair renders a solution to the 
constraints on Yii\B{pi) and Y,2\B{p2) but not on the cylinder X. The machinery of the 
conformal method — see section — is used then to find a solution to the constraints 
on Ei#I]2 using the pair {'^}'^,A]'j) as seeds for the procedure. Some amount of 
technical work has to be invested, once the solution {"flj^Kl^) has been obtained to 
guarantee that the solution coincides with the original ilfj^Kfj) and {ifj^Kfj) outside 
the gluing neighbourhoods B{pi) and B{p2)- 

The conditions on the initial data sets to be glued by the afore-discussed method 
seem, at first sight, relatively mild. Essentially, it is only required that the initial data 
sets are such that the domain of influence of the neighbourhoods B{pi) C Si and 
B{p2) C S2 have no Killing vectors. The latter condition can be reformulated in an 
intrinsic way in terms of the so-called Killing Initial Data (KID) — see e.g. [7^. More 
precisely, the condition for gluings is that there are no KIDs in B{pi) and B{p2) — in 
particular, the latter implies that there are no Killing vectors in the neighbourhoods. 
The condition on the non-existence of KIDs has been shown to be generic in a very 
specific sense — see [80] . 

Every step in the procedure as described in, say, references |338|ll82j is completely 
constructive. However, it would require to solve elliptic equations on manifolds with 
nontrivial topology. This complication is aggravated by the fact that, because of 
the requirement of the non-existence of Killing vectors in the initial data sets to be 
glued, it is not possible to use initial data sets with some symmetry — like axial 
one. Furthermore, the standard simplifying assumption of conformal flatness may 
complicate matters further. The reason is that the only conformally flat initial data 
with trivial topology is Minkowski data. Thus, the use of conformally flat initial data 
to glue would require gluing data which already have non-trivial topology. 

4-5. The Corvino-Schoen gluing method. 

An alternative gluing method for the construction of solutions to the constraint 
equations, which may be of relevance for numerical applications, is the so-called 
Corvino-Schoen ( CS) gluing construction. This method has been developed in 
|191[ I192[ I179[ 1180] and it allows to smoothly glue any interior region of an 
asymptotically Euclidean initial data to the asymptotic region of a slice of a stationary 
spacetime — e.g. Kerr initial data. For time symmetric asymptotically Euclidean 
solutions of the constraints, this method glues any interior region to an exterior region 
of a slice of the Schwarzschild spacetime. 

The procedure is best explained by considering the time symmetric situation in 
which one tries to glue a time symmetric solution of the constraints to an Schwarzschild 
exterior. As it is well known, for time symmetric initial data sets the constraint 
equations reduce to 



where [7^] in the last expression highlights the fact that in what follows the Ricci 
scalar will be regarded as a mapping between the spaces of 3-metrics on S and scalars 
on E. Under certain circumstances this mapping happens to be an isomorphism. That 
is, given a metric 7ij and a scalar / on S such that 



(7) 



(8) 
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then if a further scalar field g happens to be close enough to / — in a functional sense — 
then there is a metric 7^ close to "fij such that 

^'Mli,]=9- (9) 

The CS method makes use of precisely this property of the Ricci scalar mapping. 

Let "fij be a metric defined on a asymptotically Euclidean manifold 3-manifold S 
satisfying the time symmetric constraints. Consider a subset /C C S which is obtained 
by removing from S one of the asymptotic ends. The boundary of dIC regarded as a 
subset of S will be required to lie in the asymptotic region and to have the topology 
of the 2-sphere. Next, introduce some radial coordinate r so that the locus of dIC is 
given by r = To, with rp a constant. It is noted that /C is not necessarily compact, as 
it could well be the case that there is another asymptotic end inside JC. Denote by A 
the asymptotic end of time symmetric Schwarzschildean data. This asymptotic end is 
characterised by four parameters. Namely the mass, and the location of the centre of 
mass, d: 

As part of the gluing construction, one connects the regions /C and A by means of an 
annular region C to obtain a new asymptotically Euclidean 3-manifold E'. A positive 
definite symmetric tensor 7^- is defined on the resulting 3-manifold by requiring it to be 
identical to 7ij on JC and to 7^^- on A. On C, it is chosen so that it interpolates smoothly 
between jij and 7^^ . Now, if '■'^•'i?[7ij] is close enough to the constant vanishing scalar 
field — this can be achieved by moving dIC, regarded as a subset of S, suitably further 
into the asymptotic region — the result about the isomorphism of the Ricci Scalar 
mapping guarantees that there is a small — again in a functional sense — tensor 6^ij 
such that 7jj + Sjij is a 3-metric and more importantly 

(3)i?[7,^- + Sj,,] = 0, (11) 

on E'. The theorems ensuring the possibility of performing the above construction 
allow to show that the support of Sjij is C, so that 7^ -|- Sjij coincides with jij on /C 

and with 7,^ on A. 

''J 

The precise details of the gluing construction require the linearisation of the 
Ricci scalar mapping. In order to obtain a differential equation to be solved, one has 
to consider the composition of this linearisation and its adjoint to render a fourth 
order elliptic problem. The non-linear problem is then solved by a Newton iteration 
method. It turns out that in a similar way to what happens in the IMP gluing, the 
presence of KIDs in the gluing region C is an obstruction to the construction. The 
situation here is actually more delicate as the mere presence of asymptotic KIDs could 
be enough to cause problems. Note that as one moves the gluing region further and 
further into the asymptotic region, the problem could worsen. In order to get around 
these complications, one uses the freedom in the choice of the Schwarzschild data 
one is gluing to. Indeed, it has been shown in [191j that by a clever choice of the 
parameters m and one can get around the problems with the asymptotic KIDs. 
Unfortunately, it seems that precisely this procedure, is non-constructive. Thus, a 
numerical implementation of the CS gluing method would have to devise an alternative 
way of determining the parameters of the Schwarzschild exterior. 

Although the above discussion has been limited to time symmetric initial data 
sets which are to be glued to a Schwarzschild exterior, it has been shown that generic 
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metrics can be glued to a wide class of exteriors — including stationary ones. In this 
case, instead of just considering the Ricci scalar mapping, one has to consider the 
general constraint mapping which sends a pair of symmetric tensors on 5 to a scalar 
and a 3- vector —see [iklllTMj . 

The gluing constructions for initial data sets are remarkable geometric and 
analytical results which have provided new insights into the structure of solutions to 
the Einstein constraint equations. Important to analyse and discuss is the hypothetical 
use of these constructions in the simulation of black hole spacetimes. There are 
various issues at stake here: of course one has the conceptual ones; but there are 
also the computational aspects. Is it possible to calculate everything one requires? As 
mentioned above, a straight-forward implementation of the CS gluing technique may 
not be possible as it contains non-constructive aspects. Nevertheless, the existence of 
the gluing constructions has, at least, an indirect connection with the discussion of the 
robustness of the simulations of black hole spacetimes. In theory, one could use the 
CS techniques on the Brill-Lindquist (BL) initial data, and construct from it a new 
initial data set which will be exactly BL inside a compact set, but in addition it will 
be exactly Schwarzschildean in the reference asymptotic end. As already mentioned 
in section 12.11 a set of conserved constants for the Schwarzschild spacetime — the 
Newman-Penrose constants — are zero. On the other hand, for BL data these are 
non-zero — see |206| . As these constants are conserved along null infinity, what one 
finds is that each spacetime will have a different radiation content close to timelike 
infinity. This could be interpreted saying that each spacetime approaches in a different 
way the asymptotic state — i.e. the Schwarzschild spacetime. 

5. Evolution formalisms 

As discussed in section 13.21 different strategies in the gauge reduction process lead 
to distinct evolution formalisms. A natural classification of evolution formalisms 
in GR is according to the type of initial value problem they solve. Following the 
discussion in section [3TT1 they can be broadly classified as being Cauchy, hyperboloidal 
or characteristic. In this review we focus our attention on the Cauchy-type formalisms 
and will briefiy touch on some aspects of hyperboloidal and characteristic schemes. 

Following reference |107j we will distinguish between evolution formalisms, 
understood as the schemes devised to deal with the reduction process — involving 
e.g. the choice of variables or the strategy of resolution of the constraints — and 
evolution systems, as the concrete PDE systems set in a particular gauge and actually 
explicitly solved during the evolution. PDE evolution systems must include the 
reduced (or main evolution) system of section [3?2l but can also incorporate additional 
PDEs — e.g. fixing the gauge system. In general terms, a line can be drawn 
between i) formalisms addressing specific mathematical issues — such as existence, 
uniqueness, completeness — and ii) schemes aiming at the numerical construction of 
explicit solutions. We focus on some evolution formalisms where interaction between 
analysts and numerical relativists has been, or is expected to be, particularly fruitful. 
Excellent reviews on different aspects of this problem can be found in the literature 
p58l 14661 [3731 [2601 [4381 [4391 [4H6l [427] . The following two sections are devoted to 
Cauchy schemes: manifestly covariant schemes in section 15.1! and to the so-called 
3 + 1 formalisms in section [5TH Non-Cauchy approaches are considered in section [?751 
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5.1. Cauchy manifestly covariant schemes: generalised harmonic (or general wave 
gauge ) formulations and ZA formalism 

A classical approach to the study of the EFE consists in casting it as a wave equation 
for which there exists a well developed theory. Using the deDonder expression for 
the Ricci tensor j219[ 1229j and keeping in the left hand side only the principal part, 
equation ([1]) is written as 

~^gf"'dpd^g^, + d^,,T,j^S^, , (12) 

where = g'"^T'^j^ and Sfj.i, does not contain second derivatives of the metric. 
Prescribing the gauge — 0, equation becomes a wave equation. Since 

□x'' = —V^ (with □ denoting the scalar wave operator associated with the metric 
this choice corresponds to the use of harmonic coordinates. These coordinates 
were early used in the first results on well-posedness |2311 1228j . More recently, the 
elucidation of the role of [245j — see also |269j — led to the introduction of general 
gauge source functions H^, leading to Generalised Harmonic (GH) or general wave 
gauge conditions 

□a;'' = -r" = -Hf" . (13) 

In the language of section 13.21 this constitutes the gauge system in this covariant 
evolution formalism, and the reduced system follows from inserting (|13p into equation 
p2p . The local Cauchy problem is well-posed for this reduced system when appropriate 
initial conditions consistent with the constraints ©-(HI) are employed — see e.g. |258| 
and references therein for details and rigour. Consistency between the reduced and 
gauge systems — namely propagation of the gauge condition = — i/^ = along 
the evolution — is controlled by the subsidiary system. This follows from Bianchi 
identities, which imply the hyperbolic equation for 

VpVgp + i?^Qp = . (14) 

Imposing the gauge condition = and the constraints ©-(jll) on the initial slice 
E — which imply dtQfi = 0, cf. |258j — the vanishing of Qp along the evolution and 
the overall consistency are guaranteed. 

This subsidiary system issue exemplifies the different nature of the problems 
an analyst and a numerical relativist must address. In the analytical setting, the 
reasoning above shows the consistency of the whole scheme. In particular one no longer 
needs to consider the gauge system while studying the properties of the reduced 
system. The situation gets more complicated in the numerical setting, since numerical 
errors trigger constraint violation modes that invalidate the analytic argument — see 
also |260[ I261| for a discussion of the consequences of the intrinsic non-linearities in 
equation (HI]). The Z4 formahsm proposed by Bona et al. [lOTl flOSl fT09l [TTT] is 
particularly interesting in this context. Developed as a full spacetime generalisation 
of previous numerically motivated 3-1-1 schemes |1101 1112] 1106) . its manifest covariant 
formulation permits a general analysis in the lines of the harmonic schemes. The 
key elements in this formalism are the introduction of a spacetime form as a new 
dynamical field, and the modification of the EFE to 

i?p. + VpZ, + v,Zp = 87r(^r^,- ir^^,^ . (15) 

The solution space of this system extends that of GR, which is recovered for = 
— in fact, the EFE is recovered under the milder (non-generic) condition of Z^j, being 
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a Killing vector. The evolution of the new field is driven by the Hamiltonian and 
momentum GR constraints ([31)- (H]). Consequently, the field can be used as a 
measurement of GR constraint violations — see section II. B. in [107j and [109j for a 
discussion involving the Bianchi identity and the constraint structure of the theory. As 
a further and key development for numerical applications, reference [299] introduces 
additional damping terms in the Z4 system — and also in the GH scheme — to drive 
the solution towards the GR solution submanifold, = 0. 

Together with the early and ongoing application in analytic studies of GR 
equations — see e.g. the work of Lindblad & Rodnianski in section 12.11 and also 
the Kerr uniqueness proof by Klainerman & lonescu discussed in section [6. 1.11 — these 
covariant formalisms have led to evolution systems with a successful application in 
numerical implementations. Harmonic coordinates have been employed in numerical 
studies of the IB VP in GR |484l I485| [53] and related evolution issues [53l [54j . or in 
the study generic singularities [269] — the latter also discusses GH coordinates. In the 
current numerical relativity scenario, it is specially relevant the evolution of binary 
black holes through the merger with gravitational wave extraction by Pretorius [424j , 
for which he employed a generalised harmonic scheme [426j — see also |427j — using 
the damping terms based on the scheme in reference [299] . First-order symmetric 
hyperbolic forms of the generalised harmonic codes have been developed in [380j — cf. 
also [2281 17] — where it is studied the appropriate choice of gauge source functions 
for preserving hyperbolicity of the whole evolution system. The latter includes now 
not only the reduced system, but also the H'^ evolution equations. 

5.2. Cauchy 3+1 formalisms 

The Cauchy IVP involves the use of 3+1 spacetime spatial slicings. This applies 
in particular to the case of the harmonic systems in the previous section. We now 
consider a family of evolution formalisms developed in the 3 + 1 framework of GR 
|2141 13761 I231|, [35l 1509] . which provides the natural setting for the canonical or 
Hamiltonian approach to GR — notably inspired by GR quantisation considerations 
[89] I218| [35] 12161 [37] 1488] — and is very extended in its application to numerical 
relativity —e.g. [STUl [iSHl [5151 HHZ]- In this 3 + 1 framework, in addition to the 
splitting of spacetime in terms of a foliation by spatial hypersurfaces (St), the basic 
variables are fields naturally living on those slices Ef. In this section, we use the 
expression "3 + 1" in this restricted sense, customary in numerical relativity. We 
denote the lapse function, measuring the proper distance between slices in the 
normal n'' direction, by N and the shift vector, governing the choice of coordinate 
system on each slice, by The evolution vector is then written as t^^ = Nn^^ + (3^ . 
In this 3 + 1 decomposition, the EFE splits in two sets. The first one consists of the 
constraints in equations ©-([l]) whereas the second one — together with the kinematic 
relationship between the 3-metric and the extrinsic curvature Kij — is frequently 
referred to as the ADM evolution system (but is actually due to York |518j ): 

{dt - C(}) 7y- = ^2NK,, , 
{dt - Cf}) = -D.DjN 

+ N i^^^R,, + KK,, - 2if + 4^ [{S - E)-i^j - 2S^,]) . (16) 

The ADM system is not uniquely defined, since terms homogeneous in the constraints 
([3|)-(|4]) can be added. Moreover, it becomes an actual PDE evolution system once 
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prescriptions for the lapse N and shift /3* are actuaUy supphed. The resulting evolution 
systems have suffered from the presence of numerical instabilities and the closely 
related violation of the constraints. A variety of 3+1 formalisms have been developed 
along the years with the aim of coping with these problems — see |466| for a (non- 
updated) review. Depending on the strategy adopted to handle the constraints, 
evolution formalisms can be broadly classified as: i) free evolution formalisms if 
constraints ([3])-(l4|) are not enforced during the evolution, and ii) constrained evolution 
formalisms, that can be partially- or fully- constrained, according to the prescription 
of some or all the constraints, respectively. 

The underlying intrinsic hyperbolicity |260j of the EFE is an ingredient exploited 
in the task of obtaining a specific PDE evolution system. However, this does not 
mean that all the resulting evolution systems are necessarily hyperbolic. A clear 
example is provided by constrained schemes that, due to the elliptic nature of the 
constraints, lead to mixed elliptic-hyperbolic systems — see however |274j . Among the 
many possible criteria, in the present general discussion we will classify 3-1-1 evolution 
systems according to the presence or absence of elliptic equations in the system. 

5.2.1. Evolution systems not containing elliptic equations. A significant effort has 
been devoted to casting the EFE — under the appropriate gauge conditions — as a PDE 
hyperbolic system. The reason behind is the analytic control these systems provide on 
existence, uniqueness and stability issues [4381 1439] — see LeFloch's IHP talk for some 
examples illustrating the richness of this approach in the interplay between non-linear 
PDEs and geometry. This is the case of strongly hyperbolic systems — with their well- 
posedness results on the IVP — and of symmetric-hyperbolic systems, where energy 
estimates enhance the analytic control particularly in the IBVP permitting the use of 
maximally dissipative boundary conditions. In the generalised harmonic formalisms 
discussed in the previous section [5TT] hvperbolicitv is manifest. With regards to systems 
derived from the ADM formulation 466J, hyperbolicity analysis is more complex due 
to the loss of covariance. This also applies to certain 3-1-1 formalisms, that can be 
derived from covariant formalisms by performing a 3-1-1 splitting of the fields — e.g. 
3-1-1 versions of the ZA scheme in [107j . A first class of systems in the 3-1-1 context 
consists of first order — in time and space — evolution systems. Among these, the KST 
system f358' provides a very general class from which, for example, the Fritelli-Reula 
p63, 264^ or the Einstein-Christoffel [8] systems can be recovered. The Bona-Masso 
systems jllOl 1112] include some additional variables — expressible in terms of the 
connection coefficients of the three metric, as the BSSN system below — and can be 
related by variable reduction to KST systems |106j . A second class of 3-1-1 systems 
has been studied in which the second-order-in-space aspect of the ADM system is 
explicitly maintained. In this sense, reference [297] presents a very general second- 
order evolution system that includes the successful BSSN [4641 [75] or the NOR [40 3 j 
systems. In this context of second-order systems, we also highlight the work in [106j . 
as a predecessor of the ZA formalism discussed in section [5Tl We refer the reader to 
reference [466] for the presentations of a third subclass, the asymptotically constrained 
systems aiming at the construction of an evolution system where the constraint surface 
is an attractor — cf. in particular the discussion of the \- system and adjusted system 
approaches. Many evolution systems can be obtained from the previous systems — 
or closely related ones — by adopting distinct prescriptions for the lapse and shift. 
An exhaustive account of all the resulting 3-1-1 systems is beyond the scope of this 
article. Assessment of analytic well-posedness — strongly dependent of gauge choices — 
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and of the relations between first and second-order formulations, can be found in 

[Ml una USD IM]. 

5.2.2. Mixed elliptic-hyperbolic systems. Constrained formalisms are not the only 
avenue to mixed elliptic- hyperbolic systems. There are geometric or physically well 
motivated choices for N and /3* which are elliptic in nature and thus lead to a mixed 
system. We present some of the elliptic-hyperbolic PDE systems discussed in the 
literature: 

a) Free systems. They involve elliptic gauge choices for the slicing and/or the 
shift. An example of this is provided by the first of the axisymmetric 2+1-1-1 systems 
considered in reference [441j . A paradigmatic example is presented by Andersson 
and Moncrief in where the authors establish the well-posedness of the IVP for 
a system containing elliptic equations for the lapse and the shift. These equations 
follow, respectively, from a constant mean curvature condition and a spatial harmonic 
gauge. To the best of our knowledge, no numerical counterpart has been implemented. 
Other recent analytic studies can be found in references 163, 165, 164J. 

b) Partially constrained systems. The enforcement the Hamiltonian constraint 
as an equation for a conformal factor has been widely used in axisymmetric codes 
since the eighties (64* ,476, 227^ . More recently, maximal slicing and the Hamiltonian 
constraint have been employed in [271] . whereas the third constrained 2+1+1 system 
in [441] implements maximal slicing and imposes the momentum constraint — all this, 
still in the axisymmetric setting. 

c) Fully constrained systems. Early systems of this kind were implemented in 2D 
codes j226j . both for non-rotating 463j and rotating spacetimes pj. More recently, 
axisymetric 2+1+1 codes [159, 158^ have been used in the analysis of critical collapse. 
In this line, reference [44 Ij provides again an example of constrained system of this 
type, where ill-posedness is concluded after a maximum principle analysis of the 
involved scalar elliptic equations — see discussion in section 14.3.11 Regarding the 
full 3D case, fully-constrained schemes have been discussed in [1131 [9] and Moncrief 's 
IHP talk [1]. The fully-constrained scheme proposed by the Meudon group |113j 
makes use of a 3+1 conformal decomposition of spacetime based on slices reaching 
spatial infinity i". The elliptic part of the PDE system includes equations for the 
lapse, the conformal factor and the shift, following from maximal slicing if = 0, 
the Hamiltonian constraint and a combination of the momentum constraint and a 
Dirac-like gauge, respectively. Only two hyperbolic scalar modes are evolved, and 
the rest of the components are reconstructed by using the Dirac gauge — cf. Novak's 
IHP talk [T] — and a unimodular condition on the conformal metric. The scheme 
has been numerically tested for gravitational wave spacetimes. Regarding the scheme 
proposed in Moncrief 's IHP talk [1], it is ultimately devised to bypass the limitations 
in the extraction of gravitational radiation at finite distances by using instead the 
natural boundary of the problem — i.e. future null infinity. Even though important 
features are shared with [113| — use of a conformal 3+1 decomposition, resolution of 
Hamiltonian and momentum constraints along the evolution, elliptic gauges for the 
shift — it also represents a sound shift. Namely, it is not a Cauchy formalism but 
rather of hyperboloidal type: hyperboloidal 3-slices reaching up to J^"*" are chosen by 
means of a constant mean curvature K condition, implemented as an elliptic equation 
for the lapse. Finally, it is worth noting an interesting proposal within this scheme 
for the determination of the conformal class representative j^i, by solving an elliptic 
Yamabe equation. 
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5.2.3. Current status of 3+1 systems. Most numerical groups make use of codes 
based on 3 + 1 formalisms derived from the ADM scheme. Codes based on harmonic 
formulations have produced excellent results |424l I483j , but their use is still limited 
to a smaller part of the community. Among the 3+1 formalisms, the mainstream 
is represented by free systems not involving the resolution of elliptic equations, and 
they have provided the longest lived evolutions — in particular BSSN. In comparison, 
mixed elliptic-hyperbolic systems have offered limited applicability or are still at a 
preliminary stage — cf. |113j and Moncrief et al. scheme [1 . On the one hand, 
constrained systems are in principle expected to provide a better control on instabilities 
related to violation of the constraints. On the other hand, well-posedness in mixed 
systems is difficult to establish and, in particular, characteristic fields in the hyperbolic 
part are difficult to determine since part of the dynamics is encoded in fields solved 
through the (non-causal) elliptic part. In addition, non-uniqueness issues in the XCTS 
system discussed in section l¥. 3 . 1 1 could have strong implications on the well-posedness 
of some fully-constrained evolutions schemes enforcing a maximal condition for the 
slicing. The numerical consequences of this non-uniqueness in the elliptic sector are 
still unclear and behaviour near the critical value of the parabolic branching discussed 
in section 14.3.1 1 can be very dependent on the details of the numerical implementation. 
As a significant example, authors in reference j444j have renounced to the use of mixed 
systems in [441] . opting for a 2+1+1 version of the ZA formalism. The lesson here 
contained should certainly not be underestimated. At the same time, it should not 
preclude the development of such an alternative line of research that can only result 
in overall benefits. 

5.3. Other approaches to evolution 

We proceed to discuss some alternative approaches to the Cauchy formulation. We 
focus our attention on some aspects of characteristic evolutions and later on a 
particular version of the hyperboloidal problem. 

5.3.1. Evolution formalisms based on a Characteristic approach to GR. The 
characteristic initial value problem presented in section 13.11 provides an avenue to the 
numerical construction of spacetimes — alternative to the Cauchy approach — that 
is particularly well suited for the study of gravitational radiation. Since the pioneer 
work by Bondi and Penrose |115[ I116[ I412j on the characterisation of gravitational 
radiation in terms of null hypersurfaces, developments based on the use of Bondi- 
Sachs metric |116| I447j have resulted in successful numerical evolution — e.g. [96] — 
and radiation extraction — e.g. |55j — in a number of non-trivial cases, including 
single black hole spacetimes. An excellent account can be found in the article 
review by Winicour 515j. The main drawback of this approach is the need to 
cope with caustics formation during the evolution. A second obstacle lies in the 
difficulty of encoding physics. Alternatively, mixed Cauchy-characteristic formulations 
[94] (Cauchy-characteristic matching) offer a compromise between stable evolutions 
and accurate treatment of gravitational radiation. Full application to astrophysical 
scenarios seems currently limited as compared to the standard Cauchy approach. On 
the other hand, regarding applications in more geometric settings — in particular, the 
study of global solutions — characteristic or mixed Cauchy-characteristic approaches 
provide a stimulating framework for the fruitful collaboration between geometers and 
numerical relativists. In this sense, we highlight the ongoing activity in this line of 
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research — see e.g. [HSl I434[ I285| and references therein. 

5.3.2. Conformal field equations. [§| The conformally regular approach in numerical 
relativity is based on analytic studies by Friedrich ^242. .2441 1247j and started with 
numerical studies by Hiibner and Frauendiener [MS IMZl IMi 1131 HSl 
12341 1235] . In this approach one solves numerically a hyperboloidal initial value 
problem for the conformally regular field equations — for reviews see (235. ,236, 332] . 
The use of hyperboloidal foliations is promising as hyperboloidal surfaces combine 
advantages of Cauchy and characteristic surfaces. Instead of approaching spatial 
infinity as Cauchy surfaces do, they reach null infinity which makes them suitable 
for radiation extraction. Contrary to characteristic surfaces, these spacelike surfaces 
are as flexible as Cauchy surfaces and can be used in numerical calculations within 
the 3+1 approach. A difficulty with the conformally regular field equations is that 
equations are significantly larger than usual formulations of Einstein equations. Due 
to the large number of constraint equations, numerical errors require a stronger control 
on constraint propagation properties of the system. As there is not enough numerical 
experience with the equations, one cannot use established methods easily to deal with 
the encountered instabilities. Another difficulty is that the equations include, among 
others, evolution equations for the conformal factor. The representation of null infinity 
depends, in general, on the solution, so that the numerical boundary does not coincide 
with the conformal boundary. As a consequence, a numerical boundary treatment 
outside the physical spacetime is required and the calculation of the unphysical part 
of the conformal extension wastes computational resources. While these problems 
are not of a principal nature, they have made progress in the conformally regular 
approach difficult. The conformally regular approach as described above has been 
based on the hyperboloidal initial value problem. As such, spatial infinity is not part 
of the computational domain and one cannot calculate the maximal development of 
Cauchy data. This has various drawbacks. In the hyperboloidal approach it is not clear 
how the cut of the initial hyperboloidal surface at null infinity is related to timelike 
or spatial infinity. One would also like to be able to relate asymptotic quantities such 
as mass or momentum defined at null infinity to corresponding quantities at spatial 
infinity. An alternative conformally regular system was proposed by Friedrich that 
allows one to construct a regular finite initial value problem at spatial infinity [251j . 
The study of this problem led to analytic results concerning the applicability of the 
Penrose proposal and smoothness properties of null infinity |2551 14981 14971 15001 1501j . 
A major advantage of the underlying system, called the reduced general conformal 
field equations, is that the representation of the conformal factor is known a priori 
in terms of initial data and that the system consists mainly of ordinary differential 
equations except the Bianchi equation admitting a symmetric hyperbolic reduction. 
The finite initial value problem at spatial infinity has been studied numerically in 
|524j . The Cauchy problem could be solved for the entire Schwarzschild-Kruskal 
solution including timelike, spatial and null infinity [523j . Also, certain radiative 
spacetimes could be studied near spatial and null infinity. This approach is currently 
the only approach that allows the numerical study of both spatial and null infinity in 
a single finite picture so that one has, in principle, numerical access to the maximal 
development of Cauchy data and can calculate global spacetimes including their entire 

§ The content of his section is essentially due to A. Zenginoglu. We we are thankful for his enthusiastic 
input. 
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asymptotics. A further development of this approach might give valuable input from 
numerics to geometry as most open problems in mathematical relativity concern global 
questions. 

A valuable input from geometry to numerics is the idea of null infinity. Having null 
infinity in the computational domain would solve both the outer boundary problem 
and the radiation extraction problem in numerical relativity. Unfortunately, this idea 
could not yet be implemented in astrophysically motivated numerical calculations 
based on the Einstein equations. The main difficulty is due to the appearance of 
formally singular terms arising from conformal compactification of the metric. A 
reasonable approach is to choose a gauge for the Einstein equations in which each 
formally singular term attains a regular limit in a spacetime that admits a smooth 
conformal compactification at null infinity. The construction of such a gauge in the 
context of the characteristic approach has been known for a long time |486j . This 
method has been quite successful within the characteristic approach. The underlying 
coordinates, however, do not allow the simulation of highly dynamical strong fields 
due to formation of caustics in the light rays generating the coordinate hypersurfaces 
|515j . The implementation of this idea within the 3+1 approach has turned out to 
be exceptionally difficult numerically as well as analytically. In view of the recent 
breakthrough within numerical relativity, it seems clear that yet another formulation 
of the Einstein equations to deal with this problem cannot be regarded as a practicable 
solution. What is needed is a novel numerical treatment of null infinity that does not 
alter the successful simulation of the sources in the interior. 

An attempt to use a common reduction of the Einstein equations to study null 
infinity has been made in |402( ITO] , A scri-fixing gauge has been suggested in the 
context of a hyperboloidal initial value problem for the ADM reduction of the Einstein 
equations. In this gauge the spatial coordinate location of null infinity is independent 
of time. Such a gauge has first been constructed in [232] in the context of a frame- 
based conformally regular field equations. A similar approach is followed in Moncrief 's 
IHP talk [IJ. A scri-fixing gauge is based on a mixed elliptic-hyperbolic system. A 
constant mean curvature condition fixes the hyperboloidal foliation by an elliptic gauge 
condition for the lapse. In such attempts, an important question is how to fix the 
conformal factor. It has been shown that the representation of the conformal factor 
in terms of coordinates can be given a priori |524j . In this case, the well-posedness 
of the scri-fixing gauge has been proven. Further, each formally singular term arising 
from conformal compactification attains a regular limit at J^. Preliminary numerical 
experiments suggest that the method has some promise to be used in astrophysically 
motivated calculations. In [525j scri-fixing gauges have been studied in Minkowski and 
Schwarzschild spacetimes. A discussion of the tail behaviour of Schwarzschild is used 
to demonstrate the astrophysical relevance of the notion of null infinity for numerical 
applications. 

5.4- Some specific issues regarding evolution formalisms 

5.4- 1- Outer boundaries and the Initial Boundary Value Problem. The treatment 
of infinity poses diverse geometrical, numerical and physical challenges. More 
specifically, among the different manners of coping numerically with this issue we 
can mention compactifications of spatial infinity f424' , use of hyperboloidal slices 
together with compactification of null infinity — 326, 313 ISISI [328, 233, 232j 
I234i 1235] ■ Moncrief 's IHP talk 1 and previous section [5.3.21 — or implementation 
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of a characteristic scheme. Alternatively, one can remove infinity from the problem 
and consider the evolution of the EFE in a region bounded by an outer timelike 
boundary. This leads to the discussion of the well-posedness of the IBVP presented in 
section |31 A part of the community no longer regards at this issue as a fundamental 
problem in practise [2] , since current simulations can push outer boundaries sufficiently 
far away. However, the relevance of the topic is in particular reflected in the 
quantity of works in the subject — also manifest in the number IHP talks [1 dealing 
direct or indirectly with this topic; cf. IHP talks by Buchman, Rinne, Tiglio and 
Winicour. The work by Friedrich and Nagy [257j provided the first formalism in 
which well-posedness has been fully shown. However, its formulation in terms of 
tetrads and the Weyl tensor makes it difficult to implement using standard numerical 
techniques/infrastructures. More recently Kreiss & Winicour |368j have presented a 
system based on the harmonic formalism for which the IBVP is well-posed in the 
generalised sense, and the result has been extended to well-posedness in the classical 
sense in |367| . Generally speaking, symmetric hyperbolic systems play a critical role in 
the analysis of the well-posedness in the IBVP problem, through the use of maximally 
dissipative boundary conditions for getting rid of possible constraint violations — cf. 
[299j for an alternative approach making use of point-like damping terms. In this 
context, pseudo-differential theory of strongly well-posed systems in the generalised 
sense and techniques from semi-group theory have been used to analyse IBVP well- 
posedness — cf. gZi [MOl [13i gH [4251 UnS mi EHDI for some other recent 
references on the subject. Absorbing boundary conditions have been intensely studied 
in the last years [4091 [l36l [Ell |446] — cf- the review [450] . 

5.4-2. Gauge conditions. The discussion about numerically and physically 
appropriate choices of coordinates is a vast topic. Here, we limit ourselves to briefly 
account for some recent developments. First, regarding sources functions H^^ in GH 
formalisms, choices must be done in such a way that the hyperbolicity of the GH 
formulation is not spoiled. References [424) 14251 1426j make use of a wave evolution 
equation for the determination of H'^, but more research is needed on this specific 
topic [427]. In this sense, reference j379j has provided new gauge drivers preserving 
hyperbolicity in GH systems, and including a wider class of conditions motivated 
by successful 3-fl implementations. The further discussion of stability as well as 
the study of intrinsically GH-motivated slicings constitute an open line of research. 
In a different setting, the appropriate choice of coordinates has been crucial in the 
successful binary black holes implementations using the BSSN formalism fl4H [59] . In 
particular, modifications of the 1+log slicings and Gamma-driver condition for the 
shift, have proved to be fundamental ingredients in the moving puncture approach 
to black hole evolution, where the punctures are advected in the integration domain 
— cf. section l6.4.2l An analysis of the well-posedness of the resulting evolution system 
has been carried out in [297] . and in particular has given indications of a breakdown 
for large shifts. This binary black hole case shows the relevance of using symmetry 
seeking coordinates |272j . Recent contributions in this sense are the almost- Killing 
condition derived from a variational principle in [105i llllj , and its compatibility with 
the adoption of singularity avoiding coordinates in schemes not implementing excision 
[104]. 
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5.5. Final discussion on evolution formalisms 

Research in evolution formalisms for the Einstein equation provides a fruitful area 
for the collaboration between mathematical and numerical relativists. However, the 
analytic issues on which we have focused in this section should not shadow the full 
complexity of the problem, in particular the fact that current numerical successes are 
a direct consequence of decades of steady and systematic numerical experimentation. 
The most successful codes — e.g. generalised harmonic or BSSN — represent a 
compromise between well-posed but complicate schemes — like the Friedrich-Nagy 
system — and formalisms easy to implement but potentially ill-behaved, like the ADM 
approach. The situation could be deceiving for the analyst, but it is just a reflection 
of the overall complexity here involved. Analytic well-posedness is not a guaranty 
of numerical stability — cf. discussion in [53] — but rather a necessary condition for 
numerical convergence, once consistency of the scheme and numerical stability are 
satisfied. In sum, analytic issues are only one aspect in the problem of choosing the 
appropriate numerical evolution scheme. Together with them, the numerical relativist 
must cope with computational and physical issues |53[ . In particular, a presentation 
of the formalisms employed in numerics cannot be complete without a discussion of 
the numerical techniques involved, a task that would require a full review in itself. 

For all these reasons, the discussion about evolution formalisms is perceived by a 
good part of the numerical community as an essentially technical discussion that, in 
a good measure, has lost its critical character after the breakthroughs in the binary 
black hole problem. After decades of struggle, analytic and computational issues 
are considered to be under reasonable control and, following this line of thought, 
focus should consequently be shifted to the more physical aspects of the problem 
— cf. special New Frontiers in Numerical Relativity issue [2l and in particular the 
transcription of the final discussion. As far as astrophysical aspects in numerical 
relativity are concerned, this attitude seems to be an appropriate one — although 
surprises should be expected when addressing other physical problems. However, 
regarding the application of numerical techniques to deepen our understanding of the 
geometric structures of GR — what we have called geometric numerical relativity in 
section [1.2.21 — it seems reasonable to consider that further research in this subject is 
not only fully justified, but will probably be needed as new geometric and/or physical 
goals are set. 



6. State of the art of black holes 



If one wants to make use of numerical methods to calculate dynamical black hole 
spacetimes then, in one way or the other, one has to resort to an initial value 
formulation — see section [3l It is not obvious at all that an initial value formulation 
of the merge of two black holes is the most adequate way of discussing this problen{j]| 
— one is dealing with an idealised astrophysical problem whose initial conditions can 
not be set in a laboratory. The whole approach hinges very heavily on the assumption 
that there are certain robust aspects of the output of the numerical simulations — 
essentially the wave forms — which are not dependent on the detailed way the initial 
data are constructed. The only real dependence of the wave forms should come from 
the physical parameters of the problem — again, the masses and spin of the black 
holes — but not from things like the initial separation — as the holes rather approach 



II We thank S. Dain for discussion on this point. 
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from infinity undergoing an adiabatic process — or if tlie initial data is conformally 
flat or not. 

The issue of the robustness of certain features in the simulation of dynamical 
black hole spacetimes provides a rich arena for the interaction between Geometry and 
Numerics. This is a very challenging mathematical problem whose resolution probably 
will take long. For the time being one will have to content oneself with the evidence 
coming from the numerical simulations. Still, the assessment of this evidence is in 
itself also a challenge, and may also provide further ground for interaction. 

In defence of the use of an initial value point of view, one can say that it offers 
a framework to look at things. By framework it is being understood a repertoire of 
mathematical techniques and theorems which for example allows to state the well- 
posedness of the differential equations governing the problem, so one knows that 
the problem has been formulated in a consistent way. From a strict mathematical 
point of view, it is not known if there exist spacetimes with a black hole and 
gravitational radiation — this is essentially the problem of the nonlinear stability of 
Schwarzschild/Kerr, cf. the discussion in section [2?2l 

6.1. Mathematical black holes 

6.1.1. Rigidity and uniqueness of black holes. A cornerstone of what it has been 
called the establishment's point of view on black holes — see e.g. |415j — is the evidence 
showing that the Kerr solution describes all stationary, vacuum solutions to the EFE 
describing black holes. This uniqueness assertion is also called a no hair theorem. 
The relevance of this result stems from the fact that it characterises all possible 
asymptotic states of the general evolution of isolated systems in vacuum. Thus, for 
the community of numerical relativists it justifies the analysis of the late stages of 
black hole spacetimes by means of perturbation theory on a Kerr background — see 
e.g. [363|. 

The problem of the uniqueness of black holes has been resolved in two different 
manners. Starting from a stationary, asymptotically flat black hole spacetime it follows 
that the stationary Killing vector field must be tangent to the event horizon, £, of the 
black hole. If the Killing vector is null at the horizon — i.e. it is a Killing horizon — 
then it is possible to conclude that the spacetime is actually static [479) . For static 
black hole spacetimes, the uniqueness results state that the spacetime has to be a 
Schwarzschild spacetime — see e.g. |1751 1342( I138j . If the stationary Killing vector is 
spacelike on the horizon then from Hawking's area law one can conclude the existence 
of another vector, r;'^, tangent to the generators of the horizon which is a Killing vector 
on the horizon |310[ 13391 llSlj . This extra Killing vector field of the horizon can be 
used to define the notion of surface gravity which is of great relevance in the discussion 
of the thermodynamics of black holes. In particular if the surface gravity is non-zero, 
then the horizon is said to be non-degenerate. In the case of a non-degenerate horizon 
it is customary to assume that the horizon is a smooth null hypersurface consisting of 
two components — a bifurcate horizon — £ — U £^ intersecting at a 2-surface with 
the topology of the 2-sphere — this assumption is supported by [430j . 

Given that it is known that the Kerr solution with < a = J/M < M is 
the only stationary, axially symmetric, vacuum black hole with non-degenerate and 
connected horizon — see jl501l445] — , it is suggestive to try to extend the vector r/^ in 
the horizon to a Killing vector of the region exterior to the black hole — the domain 
of outer communication. The problem with this approach is that it requires posing 
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a boundary value problem on characteristic hypersurfaces (the horizon of the black 
hole): one constructs a wave equation for the vector t]'^ with data prescribed on £. This 
problem is ill posed — i.e. it is not possible to establish local existence and uniqueness 
by standard methods — and only admits a solution if the spacetime is taken to be 
analytic. Thus, under the hypothesis of analyticity it is possible to show that rj'^ 
extends to a Killing vector on the domain of outer communication — see [310[ 1176) . 
This type of result is known in the literature as a rigidity theorem for it shows that a 
certain assumption on the symmetry of the spacetime — stationarity — together with 
some further technical assumptions imply further symmetries — axial symmetry. The 
rigidity of black holes has been recently revisited with the particular aim of extending 
it from the 3 + 1 dimensional to the n + 1 setting — cf. Isenberg's IHP talk [1] and 
also [323] . 

The above argument for the uniqueness of the Kerr black hole has been criticised 
on the grounds that the analyticity assumption is overly restrictive: general spacetimes 
are at most smooth, and it is not obvious why a smooth spacetime will become analytic 
after it has emitted gravitational radiation. In order to get around with this problem 
an alternative approach to uniqueness has been given in |334j — see also Klainerman's 
IHP talk |T]. The idea behind this approach is that although the boundary value 
problem for wave-like equations on £ is ill posed, if one knows that solutions do 
exists, then it is possible to show uniqueness [335] . In |390[ 1391] a certain tensor — 
the Mars-Simon tensor — whose vanishing implies that the vacuum spacetime under 
consideration is locally isometric to the Kerr solution has been discussed. This tensor 
requires, for its definition the existence of a stationary Killing vector. In |334j it has 
been shown that the Mars-Simon tensor obeys a wave equation and that it vanishes 
on £ under a further technical assumption — which does not require analyticity. In 
general, one does not know whether a solution to the boundary problem for the 
wave equation satisfied by the Mars-Simon tensor with vanishing data on £ exists. 
However, if a solution exists, then it is possible to argue that the Mars-Simon tensor 
vanishes not only in a neighbourhood of the horizon but on the whole domain of outer 
communication — this is equivalent to show uniqueness of the ill posed boundary value 
problem. It is important to stress that this is a global result and not only local to 
the horizon, like the one based on the hypothesis of analyticity. Note further, that 
this proof bypasses completely the use of rigidity results, and renders the existence of 
axial symmetry as part of the main result. 

6.2. Quasi-local black holes 

In spite of the success and strength of the traditional formalisation of the black hole 
notion in the asymptotic flatness setting — cf. section [2?2l — black hole characterisation 
in terms of elements not involving global aspects of the spacetime has attracted 
significant efforts. This area of research is particularly well suited for a fruitful 
collaboration between mathematical and numerical relativists. 

6.2.1. Motivations for quasi-local Black Holes. Quasi-local characterisations of black 
holes have been approached from different communities addressing distinct problems. 
Here, we highlight geometric motivations from the mathematical community, on 
the one hand, and "practical" needs for numerical relativists, on the other hand. 
Regarding the former, the difhculties discussed in section 12.21 related to the 
understanding of the full implications of asymptotic flatness, and illustrated in the 
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absence of known examples of non-stationary "strict" black holes spacetimes, depicts 
a state of affairs which is non-satisfactory at a conceptual level. More dramatic are 
the issues raised from a numerical perspective, where the required global asymptotic 
information is simply not accesible during the evolution. Quasi-local characterisations 
of black holes are needed both for "technical" reasons — e.g. need of tracking the black 
holes along the evolution — and for physical interpretative reasons — e.g. the need to 
disentangle information about "individual black holes" in an interacting scenario. 

Other motivations have played a crucial role in the historical development of 
quasi- local notions of black holes: i) characterisation of black holes in cosmological 
settings where asymptotic flatness notions do not apply; ii) extension of black hole 
thermodynamical results to situations where global notions are not under control; 
iii) studies of black holes in the context of quantum gravity, such as the microscopic 
evaluation of black hole entropy; and iv) conceptual compliance with the consensus's 
notion of "black hole" in the astrophysical community which is completely foreign to 
global issues. In a general sense, it is fair to say that numerical and thermodynamical 
issues have been the main drivers in the development of quasi-local black holes as a 
line of research on its own. 

6.2.2. Seminal concepts and objects: trapped surfaces. As discussed in section [2.21 
the idea of a black hole as a region of no escape can be formalised in terms of the 
set points not being able of sending a signal to infinity — this, as long as we have a 
sensible notion of infinity. Local convergence of all the light rays emitted from given 
sets of points provides an alternative approach to encapture the notion of no-escape. 
This second approach does not resort to a notion of infinity and provides the rationale 
for the quasi-local approaches to black holes to be discussed in this section. 

The seminal idea behind this characterisation is the concept of trapped surface. 
This notion has played a fundamental role in the context of the singularity theorems 
|413[ 1312] . Given a surface S which is spacelike orientable and closed — i.e. compact 
and without boundary — one can unambiguously define at each point p (z S two null 
directions £^ and k^, which span the normal plane at p. That is, vector fields £^ and 
span the normal bundle T-^S. The idea behind the notion of a trapped surface S is 
that either all light rays emitted in the normal null directions do locally converge or, at 
least, light rays which would be naturally expected to expand, do actually converge. In 
the first case, the trapped character is an intrinsic property of the surface S, whereas 
in the second case some extra structure is needed to make sense of the intuition of 
naturally expanding light rays. Let us define the expansions 0^^^ and O^'^^ respectively 
associated with the null congruences £^ and k'^, as 

eW^g/^-V^C , e^'^) - g^'-V^fc, , (17) 

where is the spacelike metric induced on S from the ambient spacetime metric. 
Penrose j413j presents an intrinsic notion of trapped surface irrespective of any 
embedding. A closed surface 5 is a Trapped Surface (TS) if and only if 

g{() < , eiC^') < . (18) 

More generally, in order to cover Cosmological settings, a TS is defined as satisfying 
0W()(k) > 0, but we only discuss here the cases related to black holes. The limiting 
case in which one the expansions — say 6^^"^ — vanishes, characterises 5 as a Marginally 
Trapped Surface (MTS): 

= , 6l('=' < . (19) 
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In the context of isolated systems it is natural to consider surfaces S embedded in an 
asymptotically flat (Euclidean) spacelike hypersurface E. In this case, we can naturally 
define the outer mxW direction, say as the one "pointing" towards infinity. Following 
Hawking |3I0j we say that a surface S is an Outer Trapped Surface (OTS) if and only 
if: 

< . (20) 
The vanishing limiting case defines a Marginally Outer Trapped Surface (MOTS): 

e^^-^ = . (21) 

Dropping the global condition this last condition characterises the marginal surfaces 
defined by Hayward 315 . In order simplify the discussion, we will refer generically 
to pTjl as the MOTS condition, and the context will make it clear if global conditions 
are being taken into account or not — this is, in fact, the common practise in the 
literature, e.g. in [TS]. These are not the only attempts to classify and characterise 
the idea that a surface is trapped. See in this sense [3931 1460] for an alternative 
nomenclature, |307j for the introduction of average trapped surfaces — cf. in this sense 
also ^384j — and |46I| for a geometric formulation of the hoop conjecture in terms of 
the notion of trapped circle. 

Trapped surfaces open an avenue to control quasi-locally the notion of black hole. 
Given a spacelike hypersurface E, the trapped region 7^ is defined in 1310! as the set 
of points p G E belonging to some OTS, 5 C E. The Apparent Horizon (AH) is 
defined as the outer boundary of the trapped region Tx;. It attempts to encapture 
quasi-locally the notion of black hole horizon at the hypersurface E. A crucial result 
is the characterisation, under the appropriate regularity conditions, of the AH as 
the outermost MOTS |3I0j . Even though these quasi- local notions were originally 
developed in the mathematical studies of black holes, particularly in the context of 
the gravitational collapse and singularity theorems, they were very early employed in 
numerical implementations. 

First numerical implementations using geometric quasi-local black hole ideas. Using 
their characterisation as MOTS, AHs have played a key role in the relation between 
mathematical and numerical relativists. Firstly, for numerical relativists, condition 
(|2ip represents a tentative characterisation of "black holes" that can be efficiently 
calculated in terms of data on a compact region of a 3-dimensional slice. Important in 
this context, is the development of AH- finders, i.e. algorithms to locate AHs in spatial 
slices by solving equation (PT|) — see |490| for a review and |377) for a recent work; cf. 
also |397] for related work with a strong geometric input. Moreover, under appropriate 
energy conditions and assuming the weak Cosmic Censorship conjecture, AHs are 
geometrically defined surfaces that are guaranteed to lie inside the event horizon — see 
e.g. [310[|I86] — and therefore are causally disconnected from the rest of the spacetime. 
Although the extrapolation of this last feature from the continuum description to the 
discretised level is not straightforward — see e.g. [131j — this geometric idea has acted 
as a positive criterion in the development of black hole evolution codes — see section 
16.4.21 for a brief discussion of the different manners of exploiting this idea. 

Together with the numerical control of the black hole during the evolution, AHs 
were early applied as inner boundary conditions in the construction of initial data for 
black hole spacetimes — following an idea proposed by Unruh. In this case, condition 
(|2f p is used to complete the elliptic system defined by the constraints ©-(HI) when 
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an interior sphere has been removed — cf . section 16.4.21 for further comments on this 
excision approach. This idea were first numerically implemented by Thornburg [489j , 
later becoming a standard technique. This problem has also received attention from 
the mathematical community ,,207. 139611469] but unfortunately the interaction between 
the communities seems to be scarce. More recently, the idea of AHs inner boundaries 
has also been proposed in the evolution context by Eardley |221j . The difficulties in 
implementing these ideas exemplify the needs of interaction between geometric and 
numerical communities. See also |344[ I343j for an alternative geometric solution to 
the inner excision problem in evolutions using, instead of AHs, locally area-preserving 
evolutions of a given initial trapped surface. 

The need of a spacetime point of view. In spite of their conceptual and practical 
interest, AHs do not provide a spacetime characterisation of black holes. The AH 
notion depends on the spacelike slice E and this feature limits its applicability. The 
latter has been illustrated in examples given by Wald & Iyer [506' of Schwarzschild 
slicings where Cauchy hypersurfaces come arbitrarily close to the singularity but do 
not contain any OTS. In particular, if we understand the evolution of an AH, say 
iSo, as the "hypersurface" formed by piling up the AHs — denoted by St — found 
in the 3-slices Ei used in a Cauchy evolution — that is, D St — the resulting 
AH-worldtube IJ^ St explicitly depends on the chosen 3+1 slicing. Furthermore the 
geometric, dynamical and thcrmodynamical properties of this AH-worldtube are not 
under control — in particular non-continuous jumps can occur. This contrasts with 
event horizons, that are smooth null hypersurfaces under the appropriate conditions 

lani. 

An intrinsic spacetime formulation can be obtained by defining the trapped region 
without any reference to a particular slice E. In this context, the trapped region T 
of a spacetime Ai is the set of points p ^ Ai belonging to some TS, S C Ai. The 
outer boundary of this region is sometimes also referred in the literature as the AH, but 
following [313j we will denote it as the trapping boundary. This offers an intrinsic quasi- 
local notion for the horizon where no reference to asymptotic quantities is needed. A 
natural question to address consists in clarifying its relation with the event horizon 
— if the latter can be defined. As in the AH case, assuming Cosmic Censorship, the 
trapping boundary do not extend beyond the event horizon. Finding out if it actually 
extends up to it, or if it remains strictly in the interior, is an open issue which can 
be of relevance in the understanding of Cosmic Censorship Conjecture and in which 
mathematical and numerical collaboration can prove to be very useful — cf. ideas 
and results in |3131 12211 13691 14541 [57] . However, in spite of its conceptual interest 
this notion also entails "practical" difficulties. Namely, the trapping boundary is a 
complicate object to locate, in particular in a numerical evolution. And secondly, 
regularity issues may arise on this horizon |313| . 

6.2.3. Evolution of marginally trapped surfaces: guasi-local black holes. The need 
of an operational quasi-local notion of black hole horizon in an evolution has led to 
the formalisation of the idea of world-tubes of AHs or, more generally, world-tubes of 
MOTS — understood here as marginal surfaces in the sense of [313| . This is done at 
the price of loosing uniqueness. 

First steps where taken by Hajicek in [30 1[ 1302] 1303] were he introduced perfect 
horizons as non-expanding null hypersurfaces devised to model stationary black hole 
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horizons. These totally geodesic null hypersurfaces [301] were the first examples of 
world-tubes of AHs where not only the expansion 6*^^' vanishes but also the shear 
of the null congruence does. A more systematic and general approach was initiated 
by Hayward in |313| . where he introduced future outer trapping horizons (FOTH) 
to model quasi-local black holes in generic situations. A parallel line of research 
in quasi-local black holes as AH world-tubes has been developed by Ashtekar and 
collaborators. This has lead to isolated horizons (IH) [43ll47llt9l[44l [50 l [45 l l46 l l48] 
and, together with Krishnan, to dynamical horizons (DH) [40l [41]. FOTHs, on 
the one hand, and IHs together with DHs, on the other hand, offer complementary 
approaches to the quasi-local black hole problem. Whereas FOTHs have the virtue 
of providing a general single characterisation which encompasses both stationary 
and dynamical situations, the asset of IHs and DHs lies in their adaptation to the 
specific and rather different geometric structures associated with null (stationary) and 
spacelike (dynamical) AH-worldtubes. Both approaches underline different aspects of 
the problem and conceptual compatibility is guaranteed since equivalence have been 
shown in generic conditions [SU [181 1123] — see also [365] and later in this section. 
Finally, Booth & Fair hurst [1201ll23|l357j have introduced the notion of slowly evolving 
dynamical horizons to address the physical issues happening in the transition from 
the stationary to the dynamical regimes. General review articles on these quasi-local 
approaches can be found in [5H 11171 12921 1370] . We briefly review the ideas contained 
in them. 

Trapping horizons. A trapping horizon is (the closure of) a hypersurface Ti. foliated 
by MOTS [313j — this notion has applicability beyond the context of black holes, in 
Cosmological contexts. Here the MOTS condition = must be read in its local 
sense as a marginal surface. White and black hole situations are distinguished by 
the sign of 6'^'^'\ whereas the (local) outer- or innermost character is captured by the 
sign of 5k0^^^ — where 5k denotes formally a variation in the direction of the vector 
. Trapping horizons are originally formulated in a double-null foliation of spacetime 
|313| . where 5 becomes a Lie derivative. In this sense, trapping horizons are: 

i) Future if O'-^'i < 0, or past if O'-^^ > 0. 

ii) Outer if 46*'^' < 0, or inner if SkO*^^^ > 0. 

Note that the extremal case Skd^^'' = is not treated — see section l673l Inner trapping 
horizons are of interest in Cosmological scenarios. But in quasi-local black hole horizon 
settings, inner light are expected to converge — that is, 0^*^-' < — whereas outer 
light rays should converge just inside the black hole and diverge just outside — 
that is, SkO^^^ < 0. Therefore, the adequate characterisation is in terms of future 
outer trapping horizons, hypersurfaces Ti. foliated by MTSs St on which the stability 
condition (5^0*^^^ < holds. 

We comment on some fundamental features of FOTHs. Under the dominant 
energy condition, the MTS slices St of a FOTH have spherical topology — the topology 
law. Under the null energy condition: i) the vector tangent to the FOTH and normal 
to the MOTS slices can be either spacelike or null — the point- wise signature law — the 
latter case happening if both the shear cr'^^^ and T{i, tj vanish, and ii) the area of the 
MTS slices remains constant if Ti. is null and grows otherwise (area law). Property i) 
follows from SkO^^^ < and formalises the idea of AH-worldtubes as being stationary 
(null) if nothing falls into the horizon and superluminal (spatial) otherwise, whereas 
property ii) is the quasi- local counterpart of Hawking's area theorem [308. ,309, . 
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A main application of FOTHs has been the derivation of balance laws for physical 
quantities on the horizon, and its subsequent application to black hole dynamics — 
see Hayward's IHP talk fTI and below. Regarding the conceptual characterisation of 
quasi-local black holes, Hayward has also shown that under a regularity assumption 
— the existence of foliation of Ti. by limit sections [313] — the trapping boundary is 
actually a trapping horizon, in fact the outermost. 

Isolated horizons and dynamical horizons. The important geometric differences 
between null and spatial hypersurfaces suggest the use of different strategies for 
addressing specific issues in the stationary and dynamical regime. 

The framework of isolated horizons — see [S2] for a very comprehensive account — 
provides a hierarchy of geometric structures built on a null hypersurface, and devised 
in order to capture different stationarity levels of a black hole placed in an otherwise 
dynamical environment. 

a) Non- expanding horizons (NEH) provide the minimal notion of stationarity. 
They are defined as null hypersurfaces 7i of x K topology where: i) the null generator 
expansion 6''^^^ vanishes, and ii) the vector —T'^pf is future directed and causal. Using 
the (null) Raychaudhuri equation, a NEH 7i is a hypersurface foliated by MOTS and 
characterised by the vanishing of both cr'^^) = and T{£,£). The geometry of a NEH 
is given in terms of {q^jy, V), where q^^^ is the spatial metric induced on any compact 
slice S and V is a connection uniquely induced from the ambient spacetime connection 
V. The NEH condition expresses the evolution invariance of the intrinsic geometry: 
^£qfj.iy — 0, and essentially coincides with the notions of perfect horizons and null 
trapping horizons. 

b) Weakly isolated horizons (WIH) are NEHs with some additional structure 
needed for the analysis of equilibrium black hole dynamics: an equivalence class of 
null normals [P\ for which the surface gravities Ke, defined from £P'Vp£'^ = k^^^, are 
constant on Ti. — the zeroth law. A Hamiltonian analysis of the symplectic space of 
solutions of the EFE containing a WIH leads to a (Gibbs-like) quasi-local first law 
of black hole thermodynamics, providing in particular a quasi-local expression for the 
mass of a black hole with a stationary horizon. 

c) Isolated horizons (IH) are NEHs where the extrinsic geometry V is also 
invariant long the evolution: [£^,V] — 0. They represent the strongest quasi-local 
notion of stationarity and crucially allow the definition of mass and angular momentum 
— in the axisymmetric case — multipoles [IS]. The richness of the Isolated Horizon 
framework has led applications in very diverse areas of black hole physics, such as the 
microscopic evaluation of entropy in the context of Loop Quantum Gravity [38l |42| [39] , 
or the study of properties of hairy black holes and solitonic solutions in theories 
involving the coupling with additional fields — e.g. Yang-Mills, dilaton, Higgs, Proca 
and Skyrme fields [52^ . 

Regarding the dynamical case, Ashtekar & Krishnan introduced in [40l [41] the 
notion of dynamical horizons as spatial hypersurfaces Ti that can be foliated by MTSs 
— surfaces St where (fTQ]) holds. Note the difference with the FOTH definition, where 
the local outermost condition SkO'^^'^ < has been substituted by the condition on 
the spatial character of H. FOTHS and DHs are equivalent in generic circumstances 
[5Tl[ng. A FOTH for which Sid^^'^ ^ for at least one point on each MTS section, 
is a DH. For the converse, a genericity condition — namely SiO'--^^ ^ everywhere on 
Ti. — and the null energy condition are required. Examples of DHs not satisfying the 
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genericity condition, and thus failing to be FOTHs, are presented in [459] — see also 
comments in section 15751 in the context of extremal DHs. In contrast with the double 
null formalism of FOTHs, DHs are formulated in a 3+1 framework. Therefore they are 
specially suited for their use in numerical relativity. Although black hole mechanics 
relations can and have been derived in this setting, a strong emphasis is put in the 
application of the associated balance equations to the analysis and test of numerical 
simulations in the strong field regime. 

Quasi-local Black Hole thermodynamics. The extension of classical black hole 
thermodynamical results [63j from stationary spacetimes to more generic situations 
has been one of the main motivations for the development of quasi-local black holes 
— see Hayward's IHP talk [T]. Different approaches and definitions of the physical 
parameters have led to distinct existing versions of the generalised thermodynamics 
laws, whose detailed review goes beyond the scope of this article. We comment 
in general terms: i) Zeroth law: constancy of the surface gravity, defined as the 
non-affinity parameter associated with the generating null vector in the context of 
isolated horizons, it is shown to characterise a WIH [52] — and therefore holds on 
an IH. A generalised zeroth law in terms of an inequality bounding the mean of 
the {trapping) gravity is proposed in |313| . where the constancy of k is characterised 
by the saturation of the inequality — see also |407| for a general discussion on the 
surface gravity in the quasi-local context, ii) The second law, non-decrease of the 
area, follows from the spatial (or null character) of the MOTS-worldtubes, together 
with the condition 0^*^^ < 0. Therefore, the second law is built into the definition of 
quasi- local black holes, iii) The discussion of the first law, as a Gibbs-like expression 
relating the variation of the energy, area and angular momentum of the horizon, is 
more problematic due to the ambiguities in defining quasi-local physical parameters 
— cf. section 17.21 In the spherically symmetric and axisymmetric IH cases, this has 
been successfully addressed — see [50l [45] — and has lead to a first law in terms of 
variations in the space of physical states — Gibbs' version. However, the attempt of 
deriving a first law relating time variations of physical quantities (Clausius' version) 
has led to different versions according to the distinct choices of evolution vector or 
quasi-local energy. For these reasons, the status of such a version of the first law in the 
dynamical regime is unclear. In this context, the ensemble of balance or conservation 
laws obtained from the restriction of different components of the Einstein equation 
on the MOTS-worldtubes H acquire a particular relevance. We highlight the balance 
equation relating the variation of the area and angular momentum to the flux of energy 
[nmOmillSISllSIllIinillllj —see [Hg tor a discussion of the comparison between 
some of them. A set of evolution equations for the area |313[ [4T1 [52l I120[ 1293] , angular 
momentum j40l \M IHOl \I2T, '288', "STB] or the charge [STJ have also been derived. 
In particular, references [288, 292, 293 discuss a quasi-local version of the membrane 
paradigm [210|, 12111 14291 1493| based on a hydrodynamical analogy between horizons 
and viscous fluids — whose full analysis will plausibly involve the use of concepts and 
techniques from PDE theory of hyperbolic systems. 

We must emphasise that some of the ambiguities present in the fully dynamical 
regime can be controlled when considering slight deviations from equilibrium, in 
the setting of Booth & Fairhurst's slowly evolving dynamical horizons. Finally, 
regarding the impossibility of reaching extremal black holes {third law), a discussion 
of extremality in this quasi- local context has been recently developed in [1221 1118j . 
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6.2.4- Numerical implementations of dynamical trapping horizons. The implementa- 
tion of ideas from dynamical trapping horizons offers an example of a fruitful relation 
between Geometry and Numerics. The application of some prescriptions derived from 
the framework of dynamical trapping horizon have been used to extract physical infor- 
mation about numerically constructed black holes. Conversely, numerical experiments 
have provided key insights into the geometric structure of MOTS-worldtubes. 

DHs and "a posteriori" analyses. The main application of DHs in numerics is 
the physical and geometric analysis of MOTS-worltubcs located by an AH-finder 
along a numerical evolution. Krishnan's IHP presentation Ij provides an excellent 
account of this kind of analysis. Regarding the extraction of physics, DHs have 
been used to extract individual black hole masses and spins in black hole spacetimes 
[57, ,455t I146[ I371j . More specifically, deformation of black holes can be analysed by 
computing mass and angular momentum DH multipoles [455j . DHs can be employed 
to characterise the rate of approach to Kerr [455j or to study spin-orbit effects in 
binary black hole simulations [146j. Heuristic approaches can be formulated in this 
setting for determining a quasi-local linear momentum of a black hole, very relevant 
in the astrophysical study of recoil velocities after mergers [371 . From a geometric 
point of view, the results in [455j support the picture that AHs jumps occurring in 
black hole evolutions correspond to a situation in which the involved DHs are actually 
connected by a world-tube of MOTS which violates some of the FOTH conditions. 
The resulting single MOTS-worldtube presents timelike or signature-mixed sections 
where topology change is possible. This picture is also supported by analytic results 
in [119' on the Tolman-Bondi collapse — see also |406j . A prediction of this picture is 
that in binary black hole evolutions modelled by two MOTS-worldtubes Hi and Ti.2, 
right after the formation of the outer common horizon So, the latter should bifurcate 
into two MOTS-worldtubes: an exterior spacelike one Tioutor of growing area, and an 
interior one Thinner of mixed signature evolving into a timelike MOTS-worltube. This 
is actually seen in numerical experiments [455]. Under evolution, this Thinner should 
annihilate with the ghost horizons Ti.i and 7^2, but the details of this process are 
unknown. In this context, the assessment of the recent results in [483 showing the 
overlap of the inner horizons Tii and Ti.2 is of special relevance. 

DHs as "a priori" ingredient. An alternative application of quasi-local black holes 
to numerics is their use as a constitutive element in the PDE evolution system. 
This is exemplified by the use of NEHs to prescribe inner boundary conditions 
in the construction of initial data describing spacetimes containing black holes in 
instantaneous equilibrium —see e.g. [UHl EH [HIl [3481 [2081 [2921 [345]. A scheme 
for using MOTS — and not MTS — as inner boundary conditions in the context of 
constrained evolution formalisms has been presented in |347) and at Gourgoulhon 
& Jaramillo IHP talk 1 . This essentially recasts Eardley programme j221j in the 
dynamical trapping horizon framework. Its feasibility must still be assessed. 

6.2.5. Geometric Analysis. A fundamental shift in the research on quasi-local 
horizons has occurred with the application of tools of geometric analysis to the 
understanding of dynamical horizons properties. Besides the continuing efforts 
to characterise DH physical parameters |366| and derive appropriate conservation 
equations for them, some effort has been put into the use of maximum principles and 
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related notions to the study of the properties of eUiptic equations defined on horizon 
sections. 

The methods of geometric analysis has been put into work to show that, under 
appropriate conditions, FOTHs can be fully partitioned into NEH and DH sections — 
see [THJ[T23]. This result rules out the possibility of finding MOTS sections in FOTHs 
that are partially null and partially spacelike — a possibility which was not discarded in 
early works. This means that the transition from equilibrium to the dynamical regime 
happens all at once. This shows the complete equivalence, in generic circumstances, 
between the two main approaches to quasi-local black holes — at least in what regards 
equilibrium and the dynamical stages of the quasi-local horizons. 

Two results deserve special mention. Firstly, the foliation uniqueness theorem for 
dynamical horizons by Ashtekar & Galloway [51] : given a dynamical horizon 7i, its 
foliation by MTSs is unique. Secondly, a local existence theorem by Andersson, Mars 
and Simon [TS] stating that, given a 3-1-1 shcing (Et) and an initial MTS So C Eo 
satisfying a stability condition closely related to the outer condition of FOTHs, there 
exists a unique worldtube foliated by MTSs St, such that St C — at least as long 
as the stability condition is satisfied. 

Some more recent developments in the geometric study of dynamical trapping 
horizons include the refinement of previous work |18j on the stability and existence 
analysis of MOTS-worldtubes by Andersson, Mars and Simon in [T7]; the derivation 
of estimates for the curvature in MOTS with application on DHs [20]; studies by 
Bartnik & Isenberg [66] of necessary and sufficient conditions for the existence of 
DHs in spherically symmetric spacetimes; the analysis of the asymptotics of MOTS- 
worldtubes in spherically symmetric spacetimes |513j : the formulation of a conjecture 
about the peeling behaviour of DHs |461j in the context of the geometric discussion 
of the hoop conjecture; or the study of area estimates for outermost MOTS and 
the characterisation of the boundary of the trapped region [21] — see also |398[ 1148) . 
Finally, we mention an approach towards a general proof of Penrose inequality d la 
Huisken-Ilmanen [331j using a spacetime generalisation of the inverse mean curvature 
flows in terms of uniformly expanding flows [128j . These developments are very close 
in spirit to some of the possible applications of dynamical trapping horizons [T^] . 

6.2.6. General perspective. The recent numerical and geometric insights have 
enriched the research in quasi-local horizons, which was previously focused mainly on 
black hole thermodynamical aspects. From the DH existence and foliation uniqueness 
theorems [51] [18] it follows that the question about the evolution of a given initial MTS, 
So, is not well-deflned in generic situations. The initial So can evolve into different 
DHs depending on the chosen foliation — i.e. on the chosen lapse. Although this can 
be completely harmless in most "practical" numerical situations, it is conceptually 
important and was not sufficiently stressed in the numerical community. In numerical 
evolutions, black holes characterised by MOTSs world-tubes are treated as other 
standard compact objects. In particular, this assumes a well-defined unique evolution 
independent of the observer. The clarification of this point permits to turn the 
argument around, opening the possibility of setting a preferred choice of lapse function 
in terms of a geometrically single out DH — first attempts in this direction are discussed 
in [293j . In a different line of work, numerical results by Schnetter and Krishnan — see 
also |119j — have shed light on the global structure of tubes of MOTS, and on the 
changes of signature along the evolution and the geometric criteria to control them. 
Other issues, such as the study of generic dynamical trapping horizon asymptotics 
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towards the event horizon, and its possible apphcation to fundamental problems such 
as Kerr stability [12] or Penrose inequality will probably require the combination of 
analytic, geometric and numerical skills. 

6.3. Geometric inequalities involving black hole horizons 

As in the case of other physical theories, geometric inequalities in GR are often the 
reflect of a fundamental underlying physical principle. A prime example of this is the 
positivity of mass in GR ^456^ 14571 1516] . The fundamental nature of this result is 
manifest from its role in many crucial developments in GR. Furthermore, its failure 
would put the physical consistency of the theory under question. 

In this section we will briefly comment on some geometric inequalities in the 
context of black hole spacetimes which, in particular, constraint the gravitational 
collapse process. Our present understanding of the gravitational collapse is based on 
a chain of results and conjectures. First, the singularity theorems |4131 1311] 13121 [3T0] 
guarantee that the appearance of a trapped surface during the gravitational collapse 
leads to the development of a singular spacetime. Second, the singularity should be 
hidden behind a black hole event horizon so as to avoid a lack of predictability — this 
physically motivated hypothesis excluding the formation of naked singularities was 
proposed by Penrose [414] and is known as the weak Cosmic Censorship Gonjecture. 
Third, the black hole spacetime should reach a stationary state — this assumption 
is justified by the finiteness of the amount of radiation that an isolated system can 
emit. And fourth, assuming that all fields have fallen into the black hole after some 
finite time, the black hole uniqueness theorems [32 Ij then guarantee that the spacetime 
settles down to a Kerr black hole. Thus, barring some technical assumptions, if enough 
matter concentrates in a sufficiently compact region, then the system evolves to a final 
Kerr provided weak Gosmic Gensorship Gonjecture holds and a final stationary state 
is reached. 

Using a chain of heuristic arguments based on the previous standard picture — 
the so-called establishment picture — Penrose proposed 415 a lower bound for the 
total (ADM) mass of a black hole spacetime in terms of the square root of the area 
of the black hole. This Penrose inequality provides a lower bound for the black hole 
contribution to the total mass. It conjectures, in particular, a significant strengthening 
of positive mass theorems in the black hole context. In its first version — that can be 
referred to as global — this Penrose inequality conjectures that the area of any 
section of the event horizon £ satisfies Mj^-^j^^ > ^ As/Wtt. Remarkably, it can be 
formulated as a problem for initial data on a Cauchy surface S, providing a version 
local in time of the Penrose inequality. Given complete, asymptotically flat Cauchy 
data on E satisfying the dominant energy condition, the Penrose inequality — in the 
formulation of |325| — conjectures 

An.^n < IGttMI^^^ , (22) 

where Amin is the minimal area enclosing the apparent horizon — cf. [86] for an explicit 
construction illustrating the need of using the area of that minimal surface, rather than 
the outermost MOTS in E. There is a rigidity side to the conjecture. Namely, that 
the equality is only attained in the spherically symmetric (Schwarzschild) case. The 
Penrose inequality was proposed in an attempt to provide evidence of the violation 
of weak Cosmic Censorship. Other inequalities involving minimal surfaces in |276| 
— employed to provide lower bounds for areas of event horizon sections, as well as 
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upper bounds for the efficiency in the emission of gravitational radiation — were also 
constructed with the aim of providing evidence against Cosmic Censorship. However, 
growing evidence has accumulated in time on the generic validity of weak Cosmic 
Censorship [505j and the effort has shifted to the construction of a proof of Penrose 
inequality. Although the latter clearly does not imply the correctness of the standard 
gravitational collapse picture, it would actually provide a strong support for it. The 
spherically symmetric case has been proved in |388| . Huisken & Ilmanen )331j and Bray 
[127j have provided independent proofs for the so-called Riemannian case — where 
Kij = — but a general result has not yet been obtained. The Penrose inequality has 
evolved into a problem in its own right, becoming one of the important challenges in 
GR and Differential Geometry. Its alternative name as the isoperimetric inequality 
for black holes \281\ 12771 1280j underlines its intrinsic geometric relevance and might 
"lead to the importation into black hole theory of further useful ideas and techniques 
from global analysis" [177] —see e.g. [1001 IMl ISHl [Ml I3H3 13001 for references related 
to isoperimetric inequalities in this context. Reviews discussing the original Penrose 
argument, historical developments, main results and open questions can be found in 
[68l I392j . For early attempts to probe numerically the Riemannian case see |35H 1353) , 
and for some more recent numerical studies see, for example, [2091 1354] I352|, [346] . 

There exist some generalisations of the Penrose inequality which involve linear 
momentum [387] . charges inside the apparent horizon [2791 13201 I508[ — cf. also 
[389] and related works [3881 1336] — or a cosmological constant — cf. [1851 1278] for 
asymptotically anti-de Sitter spacetimes. Here we comment further on a particular 
sharpened version involving the angular momentum J in axially symmetric spacetimes 



A rigidity property — Dain's rigidity conjecture — has been explicitly formulated in 
[209] ! the equality holds if and only if the (exterior) initial data corresponds to 
(exterior) Kerr data. Defining A/ [87r {M^^^ + ^/M^^^ - ,P )] , inequahty ^ 

is expressed as < 1 and Dain's conjecture, = 1, characterises Kerr data by the 
evaluation of a single real number — see [2091 1346] for possible numerical applications. 
In case of being true, this would strengthen current results on characterisations of 
Kerr [3901 1391] and Schwarzschild [268' involving the evaluation of tensor quantities. 
The non-triviality of Dain's proposal can be appreciated in inequality (37) of 392]. 
The latter would provide the counterpart for a variational characterisation of Rcissncr- 
Nordstrom data, but has been found to be false [508j. 

A necessary condition for inequality (|23p to make sense, is the positivity of the 
quantity under the square root symbol. This leads to the consideration of mass-angular 
momentum inequalities in vacuum — in presence of matter such an inequality is easily 
violated. The data has to be subextremal. In a series of articles 201, 202, 203, 195, 204] 
Dain has proved the validity of | J| < M'^^^^ for maximal, vacuum, asymptotically flat, 
axisymmetric initial data. Equality is only reached for extremal Kerr. This theorem 
can be seen as a first step in the study of the non-linear stability of Kerr. These 
results have been discussed by Chrusciel in his IHP talk, and further developed in 
[178] . In spite of its characterisation in terms of initial data, these mass-angular 
momentum inequalities are spacetime properties. There has also been an interest to 
study some quasi-local versions of these inequalities involving, in particular, the local 
characterisation of extremality. In a first attempt, Ansorg and Petroff [STl 14191 [?7] have 
considered the substitution of the ADM mass by the Komar mass M^^^^^ evaluated 
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on the black hole horizon of stationary and axially symmetric spacetimes. This 
has led to the numerical construction of stationary configurations of a black hole 
surrounded by a matter torus where the quotient |>^|/-^^Komar could reach arbitrary 
high values [211 1419] or even M^^^^^ could become negative [37|- Instead, and in 
order to refer only to intrinsic quantities on the horizon, one could consider using 
the irreducible mass. This has led Petroff and Ansorg — in Ansorg's IHP talk [1] — 
to conjecture an inequality for axisymmetric stationary spacetimes only involving 
the area: 87r|J| < A. This conjecture has been further developed — including the 
charged case — in [32j, where use is made of the Christodoulou-Rufhni mass and it 
is analytically shown that extremality is characterised by the saturation of the area- 
angular momentum-charge inequality. In parallel. Booth & Fairhurst |122[ IllSj have 
undertaken an analysis of the local characterisation of extremality of black holes based 
on the dynamical trapping horizon framework of section [6T2l After also considering the 
use of geometric inequalities involving the area A — or alternatively, the vanishing of a 
locally defined surface gravity — they have opted for a characterisation of extremality 
in terms of the absence of trapped surfaces inside the apparent horizon. In terms of the 
owier/mner trapping horizon characterisation [313j . this reads SkO''^^ = 0, and leads to 
the introduction of a quasi-local parameter e satisfying e < 1, such that extremality 
corresponds to e 1 . DHs satisfying the genericity conditions [HI I123j referred to 
in section W% are found to be subextremal in this sense. In fact, dropping these DH 
genericity conditions is equivalent to the extremal characterization by S^O^^^ =0 — 
precisely the feature exploited in |459j to construct examples of spacetimes containing 
DHs but without trapped surfaces. 

Before concluding this section on geometric inequalities involving black holes, we 
must briefly comment on the so-called hoop conjecture. It proposes that black hole 
horizons happen whenever matter gets sufhciently compacted in all spatial directions. 
Stated in an intentionally vague manner, "black holes with horizons form when, 
and only when, a mass M gets compacted into a region whose circumference C in 
every direction satisfies C < AttM" [40 1|, I492j . In particular, the hoop conjecture 
offer an interesting example of interaction between geometry and numerics — see e.g. 
[1561 14621 11541 11551 I522j for some numerical studies. Recently, a reformulation of 
this conjecture as a genuine and mathematically sound — see also [384] — geometric 
inequality has been presented in [46 Ij — see also references therein for a review of the 
original conjecture. 

Black hole geometric inequalities offer a link between conceptual issues in GR 
such as Cosmic Censorship, positivity of mass or black hole dynamics — the latter, 
through the role of the area in the second law. Furthermore, its constraining role in 
the gravitational collapse process is of potential interest in numerical constructions. 

6.4- Binary black holes 

The binary black hole problem has been the main challenge for the numerical relativists 
in the last decades. The relevance of this problem lies, on the one hand, on its 
conceptual richness — this two-body problem provides a probe into the strong field 
non-linear regime of General Relativity — and, on the other hand, on its astrophysical 
interest — stellar and supermassive black holes constitute some of the main candidates 
for the detection of gravitational waves signals by interferometric antennae. Here, we 
emphasise the role of the binary black hole problem as a laboratory for the development 
and test of new numerical, analytical and geometrical ideas in GR. The developments 
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in the numerical evolutions of binary black hole spacetimes through the inspiral [135j , 
merger and ringdown phases together with the extraction of gravitational radiation 
|4241ll4lT[5^ have constituted a milestone. Since then, quite a large number of groups 
have succeeded in developing binary black hole evolution codes. The last couple of 
years have witnessed a rush in scientific activity. Given the volume of the associated 
literature, we refer the reader to the review [427j . We limit here ourselves to discuss 
some points of interest in the interaction between Geometry and Numerics. 

6.4- 1- Helical Killing vectors and binary black hole initial data. Gravitational 
radiation reaction drives relativistic systems into inspiral motion, circularising the 
orbits very efficiently at least for comparable mass systems |102j . Therefore, for two 
bodies sufficiently separated, it is natural to approximate the spiral orbits by closed 
circular ones. This physical image is geometrically en-captured by the existence 
of a one-parameter helical symmetry xa : ^ |240j of the spacetime whose 

infinitesimal generator will be denoted by — that is, = /i''. Near the binary 
system, this helical Killing vector is timelike. Sufficiently far away ft.^ becomes 
spacelike |240| . but a number T > exists such that the separation between a given 
point p G and its image by the isometry flow xt(p), is timelike — see also [287] for 
details. The integral lines associated with h^^ are helices. 

Helical symmetry is exact in theories with no gravitational radiation, like 
Newtonian gravity, (second order) post-Newtonian gravity and the Isenberg- Wilson- 
Mathews approximation to GR [3411 13401 1514] . Helical symmetry can be exact in 
full GR for non-axisymmetric systems with standing gravitational waves |101l I215j , 
although the spacetime cannot be asymptotically flat, in the sense that the ADM 
mass cannot be defined |282| — more precisely, it cannot have a smooth null infinity 
if there is no additional stationary Killing vector close to null infinity. Helical Killing 
vectors have been used in the numerical literature to model slow-motion adiabatic 
configurations of binary systems |5121I5111 15101 fi^ . This leads to the study of mixed- 
type PDEs. The light cylinder, characterised by the null character of the helical vector 
, separates an inner domain where the PDE is elliptic from an outer one where it 
becomes hyperbolic — cf. also the so-called periodic standing-wave approximation 
[231 [I30l [IHI Eia [77] and other recent numerical works [2391 [521]. Further theoretical 
developments making use of the helical Killing symmetry are given in |362j and |81j . 
In particular, the work by Klein — cf. his IHP talk yp — aims at setting a consistent 
framework for numerics by taking full advantage of the presence of a Killing vector. 
He studies the EFE in the presence of a helical Killing vector for a vacuum spacetime 
with two disconnected Killing horizons of spherical topology. The use of a projection 
formalism on the space of orbits of the Killing vector, permits then to cast the problem 
in terms of the 3-dimensional gravity equations with a 5i(2,]R)/50(l, 1) sigma model 
as material source. 

We must briefly comment on current constructions of binary black hole initial 
data. In the context of the XCTS construction — see section 14.31 — the choice of the 
evolution vector as an (approximate) helical Killing vector, i.e. = /i^, has been 
used in the literature for motivating the Ansatz 7^^^ ~ Q ~ K for the free XCTS 
data. The helical Killing vector idea was applied for the first time to the construction 
of binary black hole initial data in |2901 1295] and has subsequently been used in 
|1881 [28l [29l 1151j . All these data make use of an excision technique to deal with 
the singularity — see section [6.4.21 Problems arise when combining XCTS and helical 
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symmetry with a punctme approach to initial data, as shown in [304] . For this reason, 
punctured initial data make use of the extrinsic curvature approach — see section — 
which turns out to be difRcuh to reconcile with a helical symmetry idea. However, 
due probably to their relatively large initial separations, linear superpositions of two 
Bowen-York solutions — firstly used in [73, 58, 30] — have led to good results in recent 
punctured binary black evolutions in [Mji62ji503] and [I4 H [T42 l [T43 l fT44] . 

6.^.2. On the issue of moving punctures. A key issue in numerical black hole 
simulations is the manner in which black holes are modelled in the calculations. The 
two main techniques employed arc: excision, where a spatial neighbourhood of the 
singularity is removed from the numerical grid on the initial hypersurface and then 
subsequently, also on each spacclikc hypersurface constructed during the evolution 
[489, .458. 26 , .424, .475, ,453, .483 ^ —see also 385, .386J for particular analyses of global 
existence in this setting; and punctures, where one begins with punctured initial data 
in the sense discussed in section 14.11 and then evolves the data including the singular 
point —see for example [IM ISL llE iZl [Mi [113 Ell SSI IM] • A third ahernative 
replaces the singular black hole interior at each spacelike slice by a regular one, leading 
to the idea of stuffed black holes [331 [31] that has been recently brought back — see 
|1311 1225] . All these approaches — excision, punctures and stuffing — rely on the 
intuitive idea that no information escapes from the black hole interior, and thus, 
assume from the onset some sort of Cosmic Censorship. In particular, a detailed 
analysis of this last — non-trivial — point has been undertaken in |131j . 

There has been a recent interest in understanding some geometric aspects of the 
moving punctures picture — see e.g. |305| I306| I49H [l33] . However there seems to be 
no fully consensed view on why the puncture method works, and what happens with 
the punctures during the evolution. In particular, it is of interest to know if they still 
represent a compactified infinity. More importantly, it is crucial to know whether the 
method relies on numerical errors near an under-resolved puncture, in which case it 
could fail when probed at higher resolutions or if the evolutions are let to run much 
longer than up to now. 

Assuming the establishment point of view on black holes, the asymptotic state 
of the evolution will be described in some sense by a Kerr spacetime. So, after the 
system has evolved long enough it is natural to expect the existence of an approximate 
stationary Killing vector that could be use to drive the evolution to an eventually 
stationary slice. A crucial element in the evolution of punctured data is the use of 
the so-called symmetry-seeking choices of lapse and shift — the 1+log lapse and the 
Gamma freezing shift [5, .503. 1297j . A symmetry-seeking gauge is one that tends to 
align itself with this approximate Killing vector as the evolution proceeds. In this 
context, for a stationary slice it is understood a hypersurface such that for the chosen 
gauge-drivers, the evolution vector is parallel to the stationary Killing vector — so that 
the evolution of the spacetime effectively freezes. The discussion of stationary slices 
in the Schwarzschild spacetime has been elaborated further in [224) 14331 13051 13061 174] . 
Clearly, it is of interest to carry out an analogous analysis for the Kerr spacetime. 

The symmetry-seeking nature of the gauges used in the evolution of dynamical 
black hole spacetimes is supported by numerical evidence. According to this evidence, 
the evolution using symmetry-seeking gauges seems to pile pointwise in almost 
stationary slices of the gauge. It would be of great interest, both practical and 
theoretical to have a rigorous analysis of this symmetry seeking-behaviour of particular 
gauges. The resolution of the problem of the non-linear stability of the Kerr spacetime 
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is likely to clarify this issufl 

General results on the topology of hypersurfaces in asymptotically flat spacetimes 
— see e.g. |2651l266j — preclude the change of topology during the evolution, unless the 
slices touch the singularity — see also the discussion in [401] . section 31.6. Numerical 
evolutions — together with analytical considerations — of punctured Schwarzschild 
initial data using symmetry-seeking gauges [3051 1306j show indeed that the solution 
approaches pointwise a stationary slice that does not hit the physical singularity. 
Most importantly, the limiting slice does not reach the inner asymptotically flat spatial 
infinity, but rather ends on a cylinder of finite areal radius whose throat has an infinite 
proper distance as measured from the apparent horizon. This is interpreted as a change 
of the geometry of the puncture along the evolution — from a 1/r behaviour of the 
conformal factor near the puncture to l/v^ — that suggests the use of asymptotically 
cylindrical data to represent black holes — see also the discussion in the IHP talks by 
Hannam and O'Murchadha. These results have also been discussed in [133j . where 
it is pointed out that the puncture actually remains at the inner spacelike infinity 
along the whole evolution and, in particular, the evolution using positive lapse is 
never stationary in the limit of infinite resolution. Both apparently contradictory 
interpretations can be reconciled in terms of the differences between uniform and 
pointwise convergence: the limit to the final cylindrical stationary slice seems to be 
pointwise, but not uniform. In actual numerical simulations, the punctured region 
is under-resolved and there is no practical need to distinguish between both kinds 
of convergence. In particular, this lack of resolution is not likely to cause problems 
for finite difference codes at any reasonable — practically realistic — resolution. In 
this interpretation, the moving puncture evolution can be seen as an effective natural 
excision method [133j- It is expected that the mechanisms operating in the case of 
simulations for Schwarzschildean initial data are in essence the same ones working 
in numerical evolutions of dynamical black hole spacetimes — see comments in |133| 
about the presence of spin and momentum. 

6.4-3. Discussion. The recent intense activity in numerical evolutions of binary black 
holes has meant a very rapid advance in the understanding of the physics of binary 
black holes — with a particular emphasis on astrophysical applications. Although 
more systematic studies are needed, a large number of results are already available, 
for example concerning the effect that different configurations of physical parameters 
— such as mass ratios, spins or orbit eccentricities — have on the final state of the 
resulting black hole, in particular on the recoil velocities of the final black hole. Again 
we refer to |427| for a bibliographic account of these advances — cf. also the IHP talks 
by Lousto, Lindblom and Van Meter. 

In what concerns astrophysical applications, one should mention the complemen- 
tarity of numerical approaches with others of analytic nature. On the one hand, given 
the computational cost of the numerical simulations it is natural to explore analytic 
or semi-analytical methods to cover the whole parameter space of the problem. On 

^ It is perhaps of interest to note that in |268l a particular prescription for a lapse and a shift — 
in terms of 3-1-1 quantities — has been calculated such that if the spacetime is the Schwarzschild 
solution, then these lapse and shift coincide with those implied by the timelike Killing vector. If 
one is analysing, say, the evolution of a head-on collision of black holes where one knows that the 
asymptotic state will be the Schwarzschild solution, then it is natural to expect these lapse and shift 
— which can be evaluated in any spacetime — to be also symmetry seeking. This particular analysis 
draws ideas from the invariant characterisation of spacetimes which has been extensively developed 
in the study of exact solutions to the Einstein field equations — see e.g. | 477| . 
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the other hand, analytical approximations often offer a (quick) hint on the physics of 
the problem. The need of such a synergy between numerical and analytical efforts was 
advocated by Damour and Blanchet in their respective IHP talks [T], where different 
post-Newtonian approaches were confronted with numerical results [2131 1103j . Black 
hole perturbation theory offers a further instance of potentially fruitful combination 
of analytical and numerical approaches, for example in the analysis of the extreme 
mass ratio case — cf. Nagar's IHP talk ll. 

Astrophysical aspects of the binary black hole problem offer plenty of room for the 
collaboration between analysts and numerical relativists. In addition, the very success 
of the simulations already offer some indirect evidence of geometric issues like Cosmic 
Censorship or Kerr stability. However, the numerical control of this non-linear problem 
opens an outstanding possibility of gaining insights into some of the key geometric and 
conceptual problems of GR, if specific efforts are geared in this direction. Geometric 
ideas and results, such as Cosmic Censorship or black hole uniqueness theorems, have 
built a solid conceptual framework for the study of gravitational collapse. In turn, it 
is to be expected that strong numeric tests to some of these geometric ideas could be 
devised, offering a window to the assessment of non-linear features in GR. In current 
simulations, non-linear effects in the binary black hole have resulted unexpectedly 
mild. But, as discussed for example in Lehner's IHP talk [1], "unchartered trails" 
remain in the road and much should be learnt from studies of generic configurations 
— cf. in this sense [3811 1145] where multi-black hole encounters are studied. 

7. Miscellaneous topics 

There are some important topics in the relation between geometry and numerics 
that have not been discussed in previous sections. Here we briefly present — non- 
exhaustively — some of those aspects and lines of research. 

7.1. Critical collapse 

The study of critical phenomengfil in gravitational collapse is one of the paradigms of 
the interaction between numerical and mathematical relativity: a type of phenomena 
which was discovered by means of numerical experiments |157j which then, in turn, 
have been understood by analytical methods — see [298[ [88] for a review. In 1992 
Choptuik [157] found that it is possible to form arbitrarily small black holes in 
the process of gravitational collapse of a spherical scalar field, by fine tuning any 
parameter which affects the self-gravity of the initial configuration. Before the black 
hole is formed the spacetime becomes selfsimilar and universal, forgetting the initial 
condition except for a single length scale which determines the mass of the black 
hole to be formed. This interesting behaviour has been found in many other matter 
systems, and perturbative arguments suggest that it could be also present without 
spherical symmetry, being therefore intrinsic to the strong regime of GR dynamics. 
Finally, perfect fine tuning generates a naked singularity with infinite curvature, visible 
to observers at infinity, with important consequences for Cosmic Censorship and 
Quantum Gravity. These critical phenomena in gravitational collapse were completely 
unexpected before 1992 and are now considered the best example of how numerical 
relativity can contribute new physics to GR — see Aichelburg's and Harada's IHP 

+ The content of his section is essentially due to J.M. Martfn-Garcfa. We are thankful for his 
enthusiastic input. 
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talks for examples of ongoing research. A posteriori, they can be understood as the 
discovery of codimension 1 exact solutions in the infinite-dimensional phase space 
of General Relativity — Minkowski, black holes and stars begin global attractors of 
codimension — opening a new way of interpreting GR from the point of view of the 
theory of dynamical systems. 

7.2. Quasi-local physical parameters 

The determination of physical parameters in a finite region of spacetime is of clear 
practical importance in numerical implementations. More generally, it is of conceptual 
relevance in the theory as a whole and, more concretely, in some particular geometric 
constructions such as certain geometric flows [2771 12781 l68l I331[ I128j or in the study 
of geometric inequalities — cf. section 16.31 The determination of the amount of 
energy in a compact domain is a classical problem in GR and an exhaustive list 
of proposals regarding this goes beyond the scope of this article — see [480]. The 
relevance of this subject as a boost for the interaction between geometry and numerics 
is illustrated, for example, in the use of quasi-local dynamical trapping horizons — 
section 16.21 — for the extraction of physics in current numerical simulations. We 
highlight the attempts of defining a quasi-local angular momentum by means of 
an integral on a closed surface S involving an axial vector 0' defined on S — e.g. 
[131 iSl SHOl mi [H (mi inil IMI ISSB], see also the related work in 
and in particular the recent efforts for deriving an unambiguous prescription for 0* 
|220[ 11901 13661 1481] . Note that the divergence- free character of (p^ in not guaranteed in 
all the latter schemes, a relevant point if a slice-independent definition — i.e. intrinsic 
to S — of the angular momentum is desired. In addition to determining the magnitude 
of the spin, this vector (j)^ can be used to estimate the direction of the angular 
momentum vector — cf. [146] for preliminary results. 

7.3. Spacetime singularities and extensions of GR 

The study of the limits and extensions of GR offers stimulating perspectives of 
research. One the one hand, numerical studies of singularities [88 ] 14231 1322| 1270] — e.g. 
Cauchy horizon instabilities and mass inflation near the Cauchy horizon, naked and 
null singularities, approach to spatial singularities and analysis of BKL conjecture — 
explore the internal consistency of the theory and offer insight into new conceptual 
developments in it. Numerical implementations of higher dimensional spacetimes 
as well as the coupling of the geometry to additional fundamental fields — e.g. in 
Einstein- Yang-Mills theory — probe extensions of the theory that can be relevant in 
the quantum context. 

7.4. Some astrophysical considerations 

Ultimately motivated by the astrophysical study of relativistic compact objects, this 
article has focused on isolated systems in GR. Moreover, we have centred the discussion 
on black holes, since they offer a particularly rich context for the dialogue between 
geometers and numerical relativists. This methodological choice could shadow the 
extraordinary developments in astrophysical numerical relativity regarding systems 
with matter. Examples of the latter are the results relative to compact objects, such 
as neutron or strange stars — notably the first computations of merger of binary 
neutron stars by Shibata and UryQ [465]; see also Uryii's and Saijo's IHP talks [1] — 
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and other advances in relativistic hydrodynamics 12301 I394j . This constitutes 
presumably an extremely active area of research in next future. Moreover, the study 
of problems such as the Einstein- Vlasov system presents an analytic as well as numeric 
interest [Ml US]. 

7.5. Geometrical spacetime discretisations. 

In section [31 initial value problem formulations of GR have been presented as a 
particularly powerful approach to the construction of generic spacetimes. A radically 
different approach consists in adopting a formulation of GR which avoids the use 
of coordinates in the spirit of the Regge calculus — see e.g. [4321 1273j . For a current 
programme of research in these lines, see [23711440] — and also Frauendiener's IHP talk 
[1] — which is based on Cartan's method of frames and the use of discrete differential 
forms. 

7.6. Numerical techniques 

As we have commented in section [5] a comprehensive presentation of the evolution 
formalisms used in numerical relativity should be complemented with a discussion 
of the employed numerical techniques. Topics offering plenty of occasion for the 
collaboration between analysts and numerical relativists are, among others, adaptative 
mesh refinement, high order methods (cf. Tiglio's IHP talk) spectral methods [296] 
— including their application to evolutions [318] — and spectral elements (see also 
Maday's IHP talk), finite elements (e.g. [321 [4711 gZl STO] [36l H) finite volumes 
(e.g. [6]) or advances in high-resolution shock capturing methods. 

7. 7. Algebraic symbolic calculus 

The increase in computational capabilities and the development of powerful new 
computer algebra systems has had a great impact in many areas of GR. Many problems 
which for a long time seemed to be out reach merely on the grounds of computational 
complexities are becoming feasible. Algebraic symbolic methods were firstly used 
in GR as a tool in the general research area of exact solutions. Besides its evident 
utility in the derivation of exact solutions, computer algebra systems had a significant 
application in the metric equivalence problem — which consists in deciding whether 
two metrics given in arbitrary coordinates are locally isometric to each other. The 
equivalence of two metrics is a classical problem in Differential Geometry, and a 
solution was given by Cartan |149| . Early discussions of the equivalence problem 
in GR can be found in |126j — see also [3551 1356j . A particular implementation of 
these ideas is the system CLASSI — see |431| . Considerations involving the equivalence 
problem can be of utility in the comparison of numerical simulations. 

In recent years, algebraic symbolic calculus has been increasingly used as a 
systematic tool for analysing the analytic properties of evolution systems — e.g. in 
the construction of symmetric systems, analysing characteristics, in the construction 
of propagation systems, construction of asymptotic expansions and in perturbative 
analyses. Among the systems explicitly constructed with this purpose in mind one 
has xAct [395]. Of relevance is also the use of computer algebra systems for the 
automatic generation of computer codes like in the case of Kranc [333 . 
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8. Conclusions 

The extraordinary results in the numerical evolution of black hole binaries have had 
and will have a deep impact on the relation between mathematical and numerical 
relativists. In particular, it forces a reflection on the long term scientific objectives of 
the research of both communities. On the one hand it questions the potential pay- 
offs of certain lines of investigation, and on the other hand it offers the possibility of 
addressing problems which for technological reasons were considered out of reach. In 
particular, theory will have the unique opportunity to confront observation by means 
of the accurate simulation of relativistic astrophysical systems. Previsively, a big 
proportion of the numerical community will be involved in this endeavour. Arguably, 
the implementation of the latter will not require many geometrical insights, although 
the invariant extraction of physical content would certainly benefit from it — e.g. 
through the development of more efficient AH-finders, quasi-local characterisations 
of physical quantities, or invariant algorithms for the extraction of gravitational 
radiation. However, if one wants to study the other fundamental aspect of the theory, 
namely the structural ones — exploiting the state-of-the-art numerical possibilities — 
then evidently a geometrical perspective is fundamental. Even more, one could argue 
that the study of certain crucial geometrical aspects of GR will require numerical 
approaches to come to fruition. 

In this review we have focused on the latter geometric aspects of numerical 
relativity and also on genuinely mathematical results which, we believe, are of 
relevance for numerical applications. There is an intrinsic aesthetic appeal in the 
study of intrinsic geometric aspects of GR, but the benefit goes much further. We 
have tried to emphasise the role — and necessity — of analytical studies as guarantors 
of the internal consistency of the theory. Geometric-oriented lines research can and 
do offer general conceptual frameworks in which physically motivated problems can 
be unambiguously formulated. The gravitational collapse paradigm is an example 
of this. On a second stage, it is also undeniable that Geometry also provides 
powerful tools and insights into calculating things. Following the premise that a 
crucial aspect of "understanding a theory consists in understanding its solutions" , 
geometrical numerical relativity is one of the most powerful tools available to study 
the space of solutions of GR. Numerical GR can, potentially — i.e. assuming enough 
computer resources — solve to any desired finite precision any well formulated problem 
concerning a property of a specific solution in a concrete problem. When the problem 
involves an infinite number of solutions, then only the generic behaviour can be 
analysed. This is however where numerical GR can be most useful, finding new and 
unexpected results 0. 

We finalise by naming a number of problems for which, we believe, the interaction 
between numerical and mathematical relativists is of particular relevance — this, of 
course, notwithstanding the other issues that have been suggested in the main text. 
Cosmic Censorship remains in the eyes of most relativists the most important open 
problem in classical GR. It is to be expected that numerical investigations of global 
spacetimes will prove of utmost value in the strive towards a proof of the conjectures. 
In the closely related issue of the non-linear stability of the Kerr spacetime, numerical 
investigations have already provided information about the decays of various fields. 
This interaction is bound to become even closer in the future. A related topic, 

* We acknowledge J. M. Marti'n-Garci'a for bringing out this point. 
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the discussion of the robustness aspects — with respect to changes in the initial 
data — in the production of gravitational radiation by isolated systems will require 
strong input from numerics — even to obtain a rigorous formulation of the problem. 
The assessment of the physical relevance of certain aspects of the characterisation 
of isolated systems through the notion of asymptotic simplicity — like the peeling 
behaviour — will also require close numerical examination. Something similar can be 
said about the relevance of the construction of initial data sets by means of gluing. 
Close interaction between analytics and numerics will be required in the study of 
dynamical trapping horizons and their asymptotic behaviour close to the event horizon 
and in the relation between the global and local characterisation of the black hole 
notion — extension of the trapping boundary to the event horizon. Finally, a study of 
the solution space of GR — in the spirit of the theory of dynamical systems — using 
ideas developed in the study of critical phenomena should also be based on a close 
interaction between numerics and geometry. 
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